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We derive the quantum-mechanical master equation (gen- 
eralized optical Bloch equation) for an atom in the vicinity 
of a flat dielectric surface. This equation gives access to the 
semiclassical radiation pressure force and the atomic momen- 
tum diffusion tensor, that are expressed in terms of the vac- 
uum field correlation function (electromagnetic field suscep- 
tibility). It is demonstrated that the atomic center-of-mass 
motion provides a nonlocal probe of the electromagnetic vac- 
uum fluctuations. We show in particular that in a circularly 
polarized evanescent wave, the radiation pressure force experi- 
enced by the atoms is not colinear with the evanescent wave's 
propagation vector. In a linearly polarized evanescent wave, 
the recoil per fluorescence cycle leads to a net magnetization 
for a J g = 1/2 ground state atom. 

PACS numbers: 32.80.P - radiation pressure, 03. 75. Be - atom 
optics, 42.50.V - mechanical effects of light on atoms 



I. INTRODUCTION 

When an atom is absorbing or emitting light, its 
center-of-mass is subject to photon recoil. This phe- 
nomenon, already pointed out by Einstein in the early 
years of the century is the core ingredient of the 
atomic motion manipulation techniques that have at- 
tracted much attention during the last 20 years 
Quite paradoxically, however, this basic ingredient em- 
bodies its own limit, as it also provides crucial limit- 
ing factors to the performances of these techniques. For 
example, the minimum attainable temperatures in laser 
cooling are generally limited by the random fluctuations 
in the momenta exchanged between photons and atoms, 
that give rise to atomic momentum diffusion ; also in 
the most promising field of atom optics || , spontaneous 
emission represents a lethal threat to the coherence of 
the de Broglie waves because the atomic wave vector ac- 
quires an indeterminacy due to the random photon recoil 
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In this paper, we study the recoil effects due to spon- 
taneous emission in the vicinity of a vacuum-dielectric 
interface for an atom being reflected from an evanes- 
cent wave. This atom-optical device is one of the most 
studied realizations of a coherent mirror for atomic de 
Broglie waves §-(l| . Being based upon the interaction 
between an atom and an evanescent laser wave propagat- 
ing along the surface of a dielectric prism, it already al- 
lowed for detailed experimental investigations of fluores- 
cence rates Jl^ , p^ |, single optical pumping cycles [jl5 16 
and ground-state energy shifts Q at distances of order 
A = c/loa from the interface (u>a is the atomic tran- 
sition frequency). Spontaneous emission in the evanes- 
cent wave is also used in reflection cooling techniques 
|p^| , ^6 18 1 that have been proposed for radiative atom 
traps in the vicinity of surfaces Jl9|-p5| . From these ex- 
periments it has become apparent that a proper descrip- 
tion of the fluorescence rates and the energy levels has 
to take into account the distortion of the electromag- 
netic field due to the presence of the dielectric. A pre- 
cise theory of the atom-light interaction in the evanes- 
cent wave mirror hence touches upon the field of cavity 
QED |26[ and might even provide a model system for one 
of this field's paradigms: that the radiative properties of 
an atom are determined by the local properties of the 
electromagnetic field at the atom's position. In fact, as 
far as the 'internal' atomic dynamics (spontaneous emis- 
sion rates and frequency shifts) is concerned, this prob- 
lem has already been studied intensively, starting from 
the work of Drexhage 27 and Chance et al. |2£| and cov- 
ering a variety of geometries and materials (dielectric or 
metallic) [p9|-[l"o[. In a recent paper, Courtois et al. 
calculated the cavity QED modifications to the optical 
Bloch equations that govern the relaxation processes of a 
multilevel atom close to the vacuum-dielectric interface. 
This study showed that the relaxation due to sponta- 
neous emission is determined by the radiative damping 
rates of classical dipoles located near the interface. The 
paper was limited, however, to atoms at rest: only the 
internal dynamics was treated, whereas the 'external dy- 
namics' (recoil of the atomic center-of-mass) had still to 
be included. The present paper intends to fill this gap: 
we derive the so-called generalized optical Bloch equa- 
tions that describe both internal and external dynamics 
of an atom in the vicinity of the vacuum-dielectric inter- 
face. We actually start from a more general perspective 
and determine, for a generic cavity geometry, the mas- 
ter equation for the atomic density matrix, including the 
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center-of-mass degrees of freedom. The modifications of 
the electromagnetic vacuum field appear in this equation 
through the field correlation function, taken at two spa- 
tially separated positions. It hence turns out that the 
atomic external degrees of freedom constitute a nonlocal 
probe of the spatial correlations of the electromagnetic 
cavity field. We then follow the general procedure out- 
lined by Dalibard and Cohen- Tannoudji jfj and use the 
Wigner representation to express the spatial dependence 
of the atomic density matrix in terms of a phase-space 
quasi-distribution. In the semiclassical limit, the atomic 
Wigner function evolves according to a Fokker-Planck 
equation where appear the radiation pressure force and 
the momentum diffusion, and these quantities involve 
spatial derivatives of the field correlation function. In 
free space, this is simply a reformulation of the random 
momentum exchanges between the atom, the laser field, 
and the vacuum field p^ j. In a cavity- type geometry, on 
the other hand, photons do not rapidly escape to infin- 
ity, and the spatial structure of the field modes becomes 
important for the atomic recoil. 

We illustrate the capabilities of the Blocb-Fokker- 
Planck equation derived in this paper by focusing on 
spontaneous emission in the evanescent wave mirror. In 
particular, some unusual properties of the radiation pres- 
sure force above the dielectric interface are displayed. As 
a first example, we show that it differs from the naive es- 
timate based upon the phase gradient of the evanescent 
driving field and the spontaneous emission rate. More ex- 
plicitly, if the evanescent wave is circularly polarized, the 
radiation pressure force points into a different direction 
than the (real part of the) evanescent field's wave vector. 
This is due to the fact that the vacuum-dielectric inter- 
face partially reflects the electromagnetic field and hence 
modifies the spatial structure of the vacuum fluctuations. 
From the viewpoint of radiation reaction, the correction 
to the radiation pressure corresponds to the force exerted 
by the atom's dipole field that is backreflected from the 
interface. As a second example, we study the optical 
pumping of a J s = 1/2 atom in the vicinity of the di- 
electric. The reduced symmetry of the electromagnetic 
vacuum field implies that the average recoil per optical 
pumping cycle differs between the two Zeeman sublevels, 
even in a linearly polarized evanescent field. As a con- 
sequence, an initially unpolarized atomic ensemble splits 
into two spin components with different average momen- 
tum after a pumping cycle. For some velocity classes the 
sublevel populations then have become imbalanced, and 
the atomic ensemble shows what may be called a 'recoil- 
induced magnetization'. 

The theory outlined in this paper thus improves previ- 
ous 'heuristic' approaches to atomic recoil in evanescent 
waves that assume the atomic fluorescence to 

be distributed according to the free-space dipole radia- 
tion pattern. Our results are also relevant for radiative 
atom traps in the vicinity of material surfaces [jl9] p5[ , 
where momentum diffusion due to spontaneous emission 
may be one of the limiting factors for the temperature. 



While current atomic mirror experiments are typically 
limited to the transient regime, such traps would allow to 
study the radiation pressure force in evanescent waves in 
the long-time limit (steady state). From a more general 
perspective, the framework presented here may also be 
used to predict the center-of-mass motion of cold atoms 
in a high-finesse optical cavity with its electromagnetic 
field modes being confined both in real and frequency 
space. It is an interesting result for the domain of cav- 
ity QED that the external motion of atoms in a cavity 
provides a nonlocal probe of the cavity field correlation 
function, in opposition to internal radiative properties 
that are determined by the field correlations at the same 
point. The examples we develop demonstrate that this 
direction of cavity QED may be investigated with current 
experiments. 

The paper is organized as follows: in Sec. II, we present 
the generalized optical Bloch equations including a non- 
trivial correlation function for the electromagnetic vac- 
uum field. We focus on atoms driven by a monochromatic 
field in the low-saturation, large-detuning limit. Elimi- 
nating adiabatically the excited state, the Bloch equa- 
tions reduce to the optical pumping equation involving 
only the ground state density matrix. Passing to the 
Wigner representation, these equations take the form 
of a Fokker-Planck equation in the semiclassical limit. 
We display general expressions for the radiation pressure 
force and the momentum diffusion tensor that apply to 
any Zeeman degeneracy. The conditions of validity for 
our approach are summarized. In Sec. Ill, the example 
of the evanescent wave mirror allows us to illustrate the 
general theory. We discuss the reduced symmetry of the 
electromagnetic field correlations in the vacuum above a 
flat dielectric surface and recover, in the absence of recoil, 
the well-known fluorescence rates for this geometry. Spe- 
cializing to a scalar ground state, i.e., a J g = — ► J e = 1 
atom, we study the influence of the evanescent wave's 
polarization on the radiation pressure force and the mo- 
mentum diffusion tensor. We then consider a J g = 1/2 
ground state atom and examine optical pumping in the 
evanescent wave. The appendixes contain several techni- 
cal results that are used in the text. 



II. GENERALIZED OPTICAL BLOCH 
EQUATIONS 

The internal and external dynamics of an atom inter- 
acting with a laser field are conveniently characterized by 
a master equation for its density matrix (generalized op- 
tical Bloch equation 'G.O.B.E.'). This section is devoted 
to the derivation of such an equation in the particular 
case of a multilevel atom located close to the interface 
between vacuum and a dielectric medium. 
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A. General 



P = £00 P 



(5) 



To begin, we identify the general features of the 
G.O.B.E. at an interface. In free space, the master equa- 
tion describing the interaction of a single multilevel atom 
with a monochromatic laser field is well-known p5| , p6[ . 
Basically, its derivation proceeds in two steps. In the 
first step, one considers the evolution equation for the 
total density matrix of the system constituted by the 
atom and the electromagnetic field. In the framework of 
nonrclativistic quantum electrodynamics and in the elec- 
tric dipole approximation, this equation relies upon the 
atom-field Hamiltonian 

H = H + H R + V AL + V AR (1) 

The first term on the right-hand side of Eq.(Q) is the 
atomic Hamiltonian accounting for the internal energy 
energy of the bare atom and for its kinetic energy: 

a-^ + ^d-.-i".) (2) 

where P is the atomic momentum operator, M is the 
atomic mass, and P g and P e are the projection operators 
on the ground and excited states, respectively; the second 
term is the free Hamiltonian of the Coulomb-gauge quan- 
tized electromagnetic field; Val is the time-dependent, 
purely atomic Hamiltonian 

V AL = -D.£ L (R,t) (3) 

that describes the interaction of the atomic dipole D with 
the laser field assumed to be in a coherent state and there- 
fore described by a classical function El (r, t); and the 
last term, 

V AR = -D.E(R) (4) 

represents the coupling between the atom and the reser- 
voir associated with the vacuum quantum field E(R). 
We note that in Eq.(|l|), both fields £ L (R,t) and E(R) 
are evaluated at the location of the atom (R: atomic 
center-of-mass position operator). In the second step, 
the master equation for the atomic density matrix p is 
obtained by applying second order perturbation theory to 
the atom-reservoir interaction, and by tracing away the 
degrees of freedom associated with the reservoir. This 
yields a dynamical evolution equation where the influence 
of the reservoir is manifest through two contributions. 
The first, associated with an effective Hamiltonian, de- 
scribes the energy shifts undergone by the atomic levels 
as a result of their coupling to the vacuum field (Lamb- 
shifts). These shifts are traditionally assimilated in the 
definition of Hq, yielding the actual atomic Hamiltonian 
Ha.oo- The second contribution, p r eiax,oo, represents the 
dissipation of the atomic system due to its coupling with 
the reservoir (spontaneous emission). Finally, the free- 
space time evolution of the atomic density matrix takes 
the form 



£00 P = -T [H A .oc + Val,P] + Prelax,oo (6) 

in 

where we have introduced the free-space Liouville oper- 
ator Coa . 

We now consider an atom located in the vicinity of a 
vacuum-dielectric interface. What are the modifications 
of the master equation (^|) induced by the lower-lying 
dielectric medium? First, because of the new boundary 
conditions, the modes of El (r, t) and of the quantized 
electromagnetic vacuum field are altered and may be- 
come evanescent. It is clear that this does not affect the 
operators Hq, H r , and Val which keep the same form as 
in the free-space case. In contrast, the structure of the 
reservoir becomes modified. The contributions of Var to 
the atom dynamics (energy level shifts and spontaneous 
emission rates) are therefore expected to be different from 
the free-space situation. Moreover, as a result of the in- 
stantaneous Coulomb interaction between the atomic and 
dielectric charges, one expects a supplementary electro- 
static contribution H es to the energy level shifts. H es 
corresponds to the London- Van der Waals interaction of 
the instantaneous atomic dipole with its image in the di- 
electric medium (higher multipoles can be neglected pro- 
vided the atomic radius is much less than the distance 
between the atom and the dielectric surface). Finally, 
denoting by AHa and p re iax,int the modifications of the 
Hamiltonian and dissipative parts of the atomic density 
matrix evolution due to the interface, one obtains the 
general form of the G.O.B.E. in the presence of the di- 
electric medium 

p = Coo P + £int P (7) 

where 

£int P ~ TT" [AHa, P] + Prelax,int (8) 

in 

entirely describes the influence of the interface on the 
atomic dynamics. In particular, Ci nt p tends toward zero 
when the atom is far from the dielectric surface. The 
expression of the atomic level shifts close to a vacuum- 
dielectric interface have been presented in Ref. jll| and 
will not be discussed any further. We will therefore only 
focus on the dissipative contribution to Eq.(^). 

B. Master equation treatment of spontaneous 
emission 

We consider the relaxation processes undergone by the 
atom as a result of its coupling with the vacuum quantum 
field. As is well-known, these processes are conveniently 
described by a master equation for the atomic density 
matrix. In this section, we derive such an equation tak- 
ing into account the presence of the lower-lying dielectric 
medium. 
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1. Atom-quantum field coupling 

As stated above, the coupling between the atom and 
the quantized electromagnetic held (which is responsible 
for spontaneous emission) is described by the Hamilto- 
nian Var = — D.E(R). The atomic dipole operator D 
changes sign under parity, and therefore has only zero 
matrix elements inside the Zeeman degeneracy subspaces 
of both the ground and excited states. Furthermore, be- 
cause \g) and |e) have well-defined angular momenta, it 
is possible following the Wigner-Eckart theorem to write 
D in terms of a dimensionless, reduced dipole operator d 



T) — Vd 



(9) 



whose matrix elements contain the Clebsch-Gordan coef- 
ficients associated with the addition of the angular mo- 
menta 1 + J g — » J e . In Eq.(||), I? is a real number char- 
acterizing the electric dipole moment amplitude of the 
atomic transition. We further decompose the reduced 
dipole operator as 



d=P e dP n 



PgdP e = d^ 



(10) 



and expand d + and d = (d + )' onto the standard basis 
{u±i = T (e x ± ie y ) /y/2, u = e z } (where e x , y>z are the 
unitary vectors associated with the cartesian coordinate 
system) 



dq = d + .u q = (d q Y 



(11) 



The matrix elements of are then given by the simple 
expression 



(J e M e | \Jg Mg) = (Jg 1 Mg Q \ J £ M £ 



(12) 



where (J g \M g q\J e M e ) is the Clebsch-Gordan coef- 
ficient connecting the Zeeman sublevels \J g M g ) and 
| J e M e = M g + q). Finally, using the rotating- wave ap- 
proximation, the interaction Hamiltonian Var takes the 
more explicit form 



where 



Var = ~T>Y, « E t + d « E 'i 

q=-l 



g=-l 



(13) 



(14) 



is the positive-frequency component of the electric field 
operator. 



2. Relaxation equation for the atomic density matrix 

The total contribution p re iax = p re iax,oo + Preiax,int 01 
spontaneous emission to the time evolution of the atomic 



density matrix can be readily derived from the standard 
procedure [p9p5| |47| l outlined in Appendix |a|, where it is 
shown that in spite of the quantization of the center-of- 
mass motion, p re iax is of the familiar form, being a sum 
of two terms 

(ri| Preiax fa) = fa | { &>* (R, R) df dj , p) |r 2 ) 

+r oo C^(r 2 ,r 1 ) dj fa|p|r 2 ) d+ (15) 

where {A, B} = AB + BA denotes the anti-commutator 
between operators A and B, 

T)2 ,3 

Too = (16) 
37re n C J 

is the natural linewidth of the excited state in free space, 
and where a sum is to be taken over the i,j = x,y,z 
indices. The first line of Eq.([l5]) describes the relax- 
ation of the populations and Zeeman coherences of the 
excited state and of the optical coherences due to spon- 
taneous emission. It involves the dimensionless tensor 
C J,J fa, r 2 ), proportional to the Fourier transform of the 
electromagnetic vacuum field correlation function at the 
atomic transition frequency loq 



V 2 

rooC^fa.ra) - -j 
n 



dre jM0T {Q\Ef{Y 1 , r)^7(r 2 ,0)|0) 



(17) 



where |0) denotes the vacuum state of the field. The 
second line in Eq.(|l5[) describes the feeding of the 
ground-state Zeeman sublevels by spontaneous emission, 
yielding the expected population conservation relation 

Tr (preiax) = 0. 

It is clearly apparent in Eq.(|l|) that the effect of the 
dielectric medium on the atomic relaxation is entirely de- 
scribed by the correlation tensor C IJ (ri,r 2 ) previously 
derived by Carnaglia and Mandel (a useful represen- 
tation of C lJ (ri,r 2 ) is given in Appendix Q). Finally, 
we note that in the case where the atom is infinitely far 
from the vacuum-dielectric interface, the correlation ten- 
sor C l ' : '(ri,r 2 ) reduces to its free-space value p9[|, and 
Eq.(|l5|) transforms into its well-known expression: 

Prelax.oo fa) 

eln 

(d~.e*) e- jk ' ri fa|p|r 2 )e ik - r2 (d+.e) (18) 
where n is a unit vector and k is defined by k = (wo/c)n. 



C. Evolution of the ground state density matrix in 
the Wigner representation 

In this section, the G.O.B.E. (^) is transformed into 
a Fokker-Planck-type equation for the phase-space dis- 
tribution function of the atomic ground state. First, 
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we eliminate adiabatically the optical coherences p eg = 
P e pP g , p ge = Pg p P e and the excited state density ma- 
trix p ee = P e p P e by expressing them in terms of the 



ground-state density matrix 



PgpPg = a. This 



approximation, which holds in the limit of large laser 
frequency detunings from resonance and low saturation 
of the atomic transition, presents the advantage of re- 
ducing the atomic dynamics to a single Zeeman man- 
ifold (optical pumping equation). In a second step, 
we Wigner-transform the ground-state density operator 
a. In the new representation, a is represented by a 
(2 J g + 1) x (2 J g + 1) matrix W(r, p, t), particularly well 
suited to the investigation of the semiclassical limit of the 
atomic motion MH. 



1. Adiabatic elimination of the excited state 

In laser cooling or atom optics experiments, it is cus- 
tomary to operate in conditions of large laser frequency 
detuning from resonance and low saturation of the atomic 
transition. These conditions allow one to perform the 
adiabatic elimination of both the optical coherences p eg , 
p ge and the excited state density matrix p ee . As shown 
in Appendix [b|, this elimination amounts to the following 
substitutions. First, the dipole operators d^ 1 are replaced 
by 



d f — ► b"(R) = d [d+ • £(R)] 



£ L (r)=£ fa 



(19) 



(20) 



with b ± (R) being hermitian conjugates. The non- 
normalized, dimcnsionless vector £(r) specifies the laser 
spatial profile and polarization, while £q gives the order 
of magnitude of the electric field amplitude. 

The second replacement involves the caracteristic 
timescale of the ground-state density matrix elements: 
the spontaneous emission rate is replaced by the typ- 
ical photon scattering rate T'^ 



r i— >. r' — v 



where 



(21) 



(22) 



is the saturation parameter in the large detuning limit 
|A = cjl - lu \ > Too. 

Finally, using the usual rotating-wave approximation 
for Val an d by neglecting the influence of the atomic 
velocity on laser frequency detuning (Doppler effect), the 
G.O.B.E. (m) yields the optical pumping equation (see 
Appendix pf): 



\ri|cr|r 2 ; = — (ri| 



2M 



r 2 



-J°< ri |{C?(R),a}|r 2 ) 



+ r' 00 C«(r 2 ,r 1 )67(r 1 )<r 1 |a|r 2 )6+(r 2 ) (23) 



where 



Heff (R) = Pg A (R) Pg 

+HA' (d--r(R))(d + -«(R)) (24) 

is the effective Hamiltonian accounting for the dielectric- 
induced energy level shifts of the ground state (first term 
on the right-hand side) and the ground-state light shifts 
(last term), the order of magnitude of which is 

hA' = SAy (25) 

We also introduced in ((2^) the ground-state operator 

0(r) = CW(r,r)6+(r)&7(r) (26) 

Besides, to first order in Too/A, the connection be- 
tween a and the optical coherences and excited-state 
parts of the density matrix can be readily expressed in 
the following form (see Appendix |b]) : 

(ri|p efl |r 2 ) = (ral/vlri) 1 = 



(ri|p ee |r 2 ) = y (d + • fa)) ( ri |a|r 2 ) (d~ • ? (r 2 )) 

(28) 

where S 1 ^ is the Kronecker symbol. The validity condi- 
tions of t he expr essions given in this section are detailed 
in section II C 3. 



(27) 



2. Wigner representation of the ground state density matrix 

a. General. The Wigner representation of the den- 
sity matrix provides a particularly convenient frame- 
work for the intuitive understanding of atomic motion 
in laser light. This is because beside its intrinsic quasi- 
probability distribution character, which often enables 
to map classical pictures onto phenomena of quantum 
nature such as the atomic recoil induced by absorption 
or emission of photons, the Wigner representation is the 
best suited for taking advantage of the general charac- 
teristics of laser-cooled atomic samples to exhibit mo- 
mentum widths Ap significantly larger than the photon 
momentum fik 1 that characterizes the elementary step of 
the momentum random walk experienced by the atom 
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as a result from its momentum exchanges with the laser 
field. 

More quantitatively, the time evolution of the Wigner 
representation W(r, p, t) of the ground state atomic den- 
sity matrix 

W(r, p, t) = J d 3 s a(r: s) exp (-ip.s/ft) (29) 



c(r; s) = (r + ^s|cr|r - is) 



(30) 

can be deduced from the optical pumping equation (23). 
One thus finds after a straightforward calculation: 

dW p 
-^(r,P,*)=--^-V r W(r,p,t) + 

d 3 s 



■ip.sfh 



[ H eff(r + ±s)a{r;s) ~ a{r;s)H eff (r - ±s)] 



{G(r + is)a(r; s) + a(r; s)G(r - is)} 
4C«(r; S )b7(r+is) ( T(r; S )fe+(r-Is) 



where 



C^'(r;s) = C JJ (r-is,r + is) 



(31) 



(32) 



Because the different quantities inside the integral exhibit 
an s dependence, Eq.(pl|) clearly connects dtW{v, p, t) to 
some others W(r, p+<5p, t), which is reminiscent from the 
recoil of the atom during photon absorption or emission 
processes, hence \Sp\ » Tik. Because hk/Ap <C 1, Eq. (pi]) 
can be accurately evaluated by expanding W(r, p + <5p, t) 
up to second order in i5p [^|. A more direct way of im- 
plementing this procedure is to note that hk/Ap <C 1 
implies that the coherence length of the atomic ensem- 
ble, fr/Ap, is small compared to the optical wavelength 
A = 2n/k. This implies that the width in s of the exter- 
nal coherence function cr(r;s) is small compared to the 



scale of variation of the quantities C l J (r 



is) 



and b ± (r± is), which is of the order of A. It is therefore 
possible to expand directly the integral kernel of Eq. ( pi] ) 
up to second order in fcs before evaluating the integral. 
We now discuss more quantitatively this procedure in 
order to identify the influence of the vacuum-dielectric 
interface on the atomic dynamics. 

b. Zeroth order: internal atomic dynamics. The low- 
est (zeroth) order in the expansion of the different quan- 
tities in the integration kernel of Eq.(31) amounts to con- 
sidering the atoms as point-like particles, i.e., to treating 
the atomic translational degrees of freedom classically. 
It is therefore not surprising to end with the previously- 
established optical pumping equation pl[ for a point 
atom having a constant velocity p/M which yields, af- 
ter adiabatic elimination of the excited state and optical 
coherences: 



0[(ks)°] 



dW 
~dt 



■VrW 



M 



--{H eff {r), W}+W re i ax 



(33) 



where the effective Hamiltonian H e ff only impacts on 
the evolution of the ground state Zeeman coherences (the 
space dependence of -He// has no effect on the atomic 
motion to zeroth order in fcs ) and where 



W r . 



.Ik C ^(r;0){bt(r)bj(r),W} + 
-r' oo C^(r;0)bJ(r)Wbi(r) (34) 



accounts for departure from the ground-state through 
laser absorption (first term on the right-hand side) and 
for feeding of the ground-state by spontaneous emission 
(second term), the combined action of which yields opti- 
cal pumping. Eq. (|3~i| ) shows that the optical pumping or 
fluorescence rates are determined by the one-point cor- 
relation tensor C7*' J (r;0) = C IJ (r,r). In free space, one 
has C % J (r; 0) = S l J so the position dependence of the 
fluorescence and optical pumping rates only arises from 
the driving field profile £(r) [cf. Eq.(^9])]. On the other 
hand, close to a vacuum-dielectric interface, as will be 
shown below, one has C l '-'(r;0) = 8 t ^c % {z), so an addi- 
tional cause for a spatially varying optical pumping rates 
appears. As already pointed out in Ref. Q|, this phe- 
nomenon is directly connected to the well-known space- 
dependence of the damping rates of classical oscillating 
dipoles close to a vacuum-dielectric interface. 

c. First order: radiative and level shift-induced forces. 
To first order in fcs, where the effect of the atom-field cou- 
pling on the atomic motion enters into play, the ground- 
state G.O.B.E.'s expansion takes the form of a Liouville 
equation uncovering the force F(r) acting on the atom 



F(r) ■ V P W = 



(35) 



The force operator F(r) is actually the sum of three 
terms: 



P = jp(shift) _|_ p(dip) , p(sp) 



(36) 



where 



P (.w/t) + F (dij>)j . v p H/ = -i {V r H eff (r), V P W} 

(37) 



involves the sum of the force associated with the 
dielectric-induced energy level shifts (F<- shifi )) and of 
the dipole force (F^ dlp >) associated with the ground-state 
light-shifts (F ( - shl M and F^ are associated with the 
first and second term of Eq.([j4"|), respectively), and where 



tt( s p) - F (sp) -I- "F (sp) 

F ~ F depart ' F feed 



(38) 



G 



is the radiation pressure force, having two contributions 
arising from each term of the right-hand side of Eq. (34), 
respectively. One finds: 

F^dto ' V P W = -iftT' 00 V../4(r; 0) • V V W (39) 

where *4(r; s) is an operator involving the field correlation 
function at two points separated by s, defined as 



A(r; s)W = C 1 ' 1 (r; s) bj (r + is) W 6+(r - is) (40) 



and 



depart 



(r)-V p W = -hr t 00 \V r g(r),V p W] 



(41) 



In order to make the physical content of F^^[ more 
transparent, let us consider the simple situation where 
the driving laser field is a plane wave of wavevector q, in 
which case the operators take the simple form 



b ± (r)=/? ± e 



(42) 



where j3 Q are space-independent operators. A straight- 
forward calculation then yields 



' feed 



(r).V p Ty = ^r' oo (qC«(r;0) 



-iV s C^(r;0))f3^.V p W/3+ (43) 

Eq.(ft3|) shows that F^^d results from two contributions, 
the physical significance of which can be deduced by re- 
ferring to the well-known free space situation: the first 
term in parentheses describes the atomic recoil h<\ due to 
the absorption of the driving plane wave photons, while 
the second is associated with the atomic recoil during 
spontaneous emission. One can thus consider that the 
presence of the dielectric medium affects the quantitative 
value of Ff^j through the modification of C lJ (r; 0), but 
that it remains qualitatively analogous to the free space 
situation. 

The situation is quite different for F^^.. Considering 
Eqs. ( |4"i"| ) and (^6|), one can indeed note that whereas in 
free space, the contribution to F^^ rt only arises from the 
space dependence of the laser field, a supplementary con- 
tribution shows up in the vicinity of a vacuum-dielectric 
interface due to the space-dependence of the one-point 
correlation tensor C z,:, (r, r). In the preceding case of a 

plane wave driving field, whereas F^p-* art reduces zero in 
free space, one thus obtains a purely dielectric-induced 
contribution (cancelling in free space) of the form: 

Fffart-VpW = -M^VrC^O). [f3+^,V p W] (44) 



d. Second order: momentum diffusion tensor. To sec- 
ond order in ks, we find a Fokker-Planck equation for the 
Wigner matrix. Its diffusion term is given by 



0[(ks) 



dW 
~~dt 



d 2 W 
dpkdpi 



The momentum diffusion tensor 



D k \v)=Dl 



k,l 
epart 



D 



k,l 
Iced 



(45) 



again contains contributions from the departure and feed- 
ing terms on the right-hand side of Eq.(p4|), respectively: 



D 



k,l 

depart 



d 2 W 
dpkdpi 



h 2 T' 



d 2 g d 2 W 



16 1 dr k dn ' dpkdpi 



n k,i _ _ 

^feed — 



ft 2 !^ d 2 A 
2 ds k dsi 



(r;0) 



(46) 



(47) 



As it is well-known J50[ , this tensor results from various 
phenomena: randomness of atomic recoil processes due to 
laser absorption and stimulated/spontaneous emission, 
spatial spreading or shrinking of the atomic wavepacket 
due to the space variation of the absorption or optical 
pumping rates. Because the dielectric medium affects 
both, absorption and spontaneous emission processes, 
one expects modifications of the atomic momentum dif- 
fusion in the vicinity of the vacuum-dielectric interface. 



3. Validity conditions of the derivations 

To conclude this section, we summarize the validity 
conditions for our approach. The mose stringent condi- 
tion arises from our approximation that the excited state 
density matrix adiabatically follows the ground state den- 
sity matrix. This means first that the atoms move little 
on the scale A during the lifetime l/T^ of the excited 
state, or equivalently: 



kp kAp 



(48) 



Second, the force F e acting on the excited state must be 
sufficiently small that during the lifetime l/r^, the shift 
of the atomic momentum is negligible compared to the 
width Ap of the momentum distribution: 



Too > 



Ap 



(49) 



Under typical experimental conditions (distance z ~ 
1/fc), the force F e is at most of order |3l|] V 2 je^z^ ~ 
hkT 00 /(kz) 4: ~ hkToe. Condition j49| ) hence reduces 
to the semiclassical regime Ap ^> hk we shall suppose 
throughout this article [cf. condition (|5|) below]. 

The elimination of the optical coherences is governed 
by a different condition: the laser detuning A must be 
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larger than any other frequency scale in the G.O.B.E. (cf. 
Appendix |b|) 



|ai » r c 



VEq kp 
~' M' 



(50) 



This condition implies a small saturation parameter (|l9| ) 
and amounts to neglecting the Doppler shift of the laser 
frequency. It also allows to discard the dielectric-induced 
shift of the atomic transition frequency compared to the 
free-space detuning (see Ref. Q for more details). 

Furthermore, since oja 3> |A| (near-resonant excita- 
tion), the atomic dynamics is 'frozen' at the timescale 
of the vacuum field correlation time. This justifies the 
Markov approximation made in deriving the G.O.B.E. 
(|l5|). We note that if the conditions (^, |5(]) are re- 
laxed, one has to take into account both the ground- and 
the excited-state manifolds, with spontaneous emission 
inducing transitions between them. Thispicture is remi- 
niscent of the dressed-state description |5l| and has been 
used, e.g., in Refs. OjL9|. 



As pointed out in Subsec. [I C 2 above, the mechanical 
effects of spontaneous emission may be described simply 
in a semiclassical way if the atomic momentum distribu- 
tion varies smoothly on the scale of the photon momen- 
tum 



Ap > hk 



(51) 



Combining with condition (plq), our approach is limited 
to transitions with 3> w rcC oii = hk 2 /2M. It therefore 
fails for light atoms like He or Li, e.g., while it applies 
for heavier atoms like Na, Ne*, Ar*, Rb, Cs . . . 



III. ATOMIC MOTION AT AN EVANESCENT 
WAVE MIRROR 

In this section, we illustrate the capabilities of the ap- 
proach developed above by applying it to the motion of 
an atom in an evanescent wave mirror, in the vicinity 
of a vacuum-dielectric interface. We first examine the 
electromagnetic field correlation tensor in this geometry, 
with particular emphasis on its symmetry properties. We 
thus recover the well-known atomic damping rates above 
the dielectric interface. The radiation pressure force and 
the momentum diffusion tensor are then explicitly calcu- 
lated for a J g = — > J e = 1 (scalar) atomic transition. 
The optical pumping processes of a J g = 1/2 atom in the 
evanescent wave are also investigated. 



A. Electromagnetic field above the dielectric 

1. Vacuum field correlation tensor 

The field correlation tensor C IJ (r;s) in the vacuum 
half-space above the dielectric has been calculated by 



Carnaglia and Mandel As shown in Appendix |(], 

this tensor can be conveniently written as the sum of the 
free-space correlation tensor and an interface-dependent 
part: 



C* =C% + C£ t 



(52) 



The fre e-space correlations are given explicitly in 
Eq. ( |C19| ) , although we may deduce most of their proper- 
ties from symmetry considerations. First, due to trans- 
lational invariance, C^(r;s) = C£f(s) is independent of 
the position r and only depends on the difference vector 
s. Second, due to rotational invariance, the correlation 
tensor at the same point (i.e. s = 0) is proportional to 
the Kronecker symbol <$*'•'. Third, for s / 0, the tensor 
may be decomposed into an isotropic part proportional to 
S 1 ^ , and a quadrupolar part proportional to s's- 7 — is 2 (5 IJ . 
The coefficients of this decomposition are scalar fun ction s 
of s 2 p2[ ; we give their expansion for small s in Eq.(C19). 

Let us now apply these symmetry arguments to the 
interface-dependent part of the correlation tensor, CijL. 
We observe that the translational and rotational sym- 
metries are broken and reduce to translations parallel to 
the interface and rotations around the interface normal, 
respectively. As a consequence, we expect C^ t (ri,r 2 ) 
to depend on the distances Z\, z 2 of the interface and 
on the in-plane difference vector S|| = r|| 2 — i"||.i [the || 
subscript denotes the translational directions parallel to 
the interface (x, y components)]. More precisely, because 
the contributions to C*'- 7 arise from the partial reflection 
of the field at the vacuum-dielectric interface and the 
evanesce nt w aves present in the vicinity of the dielectric 



[cf. Eq.( |C20|) 1, C\'^ t is expected to depend only on the 



sum Z\ + Z2 of the distances from the interface. There- 
fore, we may write 



±s,r + is) 



Cmt(^; s n)- 



(53) 



As shown in Appendix this correlation tensor actually 
contains four parts having different symmetry properties: 

C^ t (z; 8,|) = c (z; sf ) + q (z; sf ) (s'W - |^) 
+fca 1 (z;s 2 )(><' i S |j - s|<S 2J ) 

+fc 2 <&(*;sj|)(«| sj - ±s 2 (<5^ - S'Wj) (54) 

The first term on the right-hand side of Eq. (|54|) provides 
a scalar contribution, whereas the second has the sym- 
metry of the Yq spherical harmonic with respect to the 
interface normal (quadrupolar part). Note that these two 
terms entirely account for the one-point correlation ten- 
sor C t w(z; 0), which is related to the modifications of the 
natural widths of the excited state Zeeman sublevels by 
the interface, as shown below. The third contribution to 
Eq. (p4[) , proportional to e IJ,fc (s|| x e z )k, displays an axial 
symmetry. Finally, the fourth term corresponds again to 
a quadrupolar part with respect to the in-plane vector Sm . 
The scalar weight functions cq, go, ax and qi are given in 
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Apppendix ^ [Eqs. (C24)] and are plotted in Fig. [| for 
s II = as a function of z. It clearly appears on this figure 
that the influence of the interface is only significant for 
distances smaller than the optical wavelength A. 
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FIG. 1. Dimensionless functions (C24) determining the 
field correlations as a function of distance z from the sur- 
face (in units of 1/fc). 

Lateral distance su = 0, refractive index no = 1.5. 

As will be shown in the following, the axial part of the 
field correlations is at the origin of nonstandard effects 
close to the dielectric surface, so a physical interpreta- 
tion of this term might be helpful. To this end, we use 
the fact that the correlation tensor C^fri,^) is pro- 
portional to the electromagnetic field susceptibility, i.e. 
the electric field Ei created at ri by a classical dipole 
located at r2, oriented along the ej axis and oscillating 
at the ato mic resonance frequency p9]J^|39| , fIo| (c/. also 
Appendix |C 3| ) . Consider now a dipole at r 2 = (0, z), os- 
cillating at the atomic resonance frequency and oriented 
perpendicular to the interface, and the field it creates at 
ri = (s||,z) (see Fig. ||). If the dipole's distance z is 
large compared to the optical wavelength, we may use 
geometrical optics to find the rays that reach the obser- 
vation point ri. As far as C\£ t is concerned, the only 
ray to consider is the one that reaches ri after one re- 
flection from the interface (the thick solid line in Fig. ||). 
The vertical orientation of the dipole implies that this 
ray must be TM-polarized. Due to the finite distance sm 
parallel to the interface, the reflected field vector at ri 
has a nonzero component parallel to the interface (actu- 
ally, parallel to sy). This construction hence illustrates 
how the reflection at the interface creates a correlation 
between lateral and perpendicular field components at 
spatially separated positions [described by the axial part 
of the correlation tensor (Q)]. If the dipole's distance z 
is not large compared to A, geometrical optics fails and 
one has to take into account a continuous distribution 
of modes that also contains evanescent waves. However, 
these modes still have the common feature of being TM- 
polarized and therefore also contribute to the axial co- 
efficient ai(z;s|) of the correlation tensor [see the Som- 
merfeld integral flC24c )1- 




FIG. 2. Illustration of the axial part of the correlation ten- 
sor: electric field created at ri by a dipole located at r-2 and 
oscillating perpendicular to the interface. 



2. Connection with the damping rates of the excited state 
Zeeman sublevels 

As shown in Ref. the atomic internal relaxation 
processes associated with spontaneous emission close to 
a vacuum-dielectric interface can be entirely described by 
means of the damping rates r^z) and T ±{z) of classical 
oscillating dipoles located at a distance z above the di- 
electric medium, polarized parallel or orthogonal to the 
interface, respectively, and such that ru ±{z — > oo) = 
Too. More precisely, one shows that the contribution 
Pee.reiax of spontaneous emission to the evolution of the 
excited state part of the internal atomic density matrix 
for an atom at rest in r reads [|41[ 

Pee,relax (^^) 

= -~ {r,| (*) d+d- + r„ (z) d+d- + T x (z) dtd-, Pee } 



By comparing Eq.(|55|) with Eq.(|15[), which yields 

r 

Pee,relax = {C* J (r, r) df dj , p ee } 



(56) 



one can readily deduce that 

= l + c (z;0)-i go (z;0) 



rx(*) 



l + c (z;0) + fgo(z;0) 



(57a) 



(57b) 



As can be checked on the expression of cq and go 
[Eqs. QC24aD and flC24b| )], these results are consistent 
with the well-known form of T\\(z) and Tj_(z) P8|]41[. 
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3. Evanescent driving field 

We focus in this paper on an evanescent driving laser 
field without polarization gradient, i.e., we consider the 
field profile 



£( r ) = Co exp (— nz+iQx) 



(58) 



where Q 2 



k 2 and £0 is a constant vector. This 



field is created by total internal reflection of a plane laser 
wave inside the dielectric (hence, Q > k). For the two el- 
ementary polarizations TE and TM of the incident wave, 
we have 



ATE) 

so 



ATM) 
so 



k 



(59) 



Note that in the TM-case, the polarization of the evanes- 
cent wave is elliptic: it approaches a linear polarization 
in the vicinity of the critical angle (k — > 0) and a circular 
one (cr~ with respect to the positive y-axis) far from the 
critical angle (k ~ Q). 



B. J a = 0^ J e 



1 atomic transition 



In this subsection we consider the simple situation of a 
scalar atom (a J g = — > J e = 1 transition) and calculate 
the radiation pressure force and the momentum diffusion 
tensor. The Wigner function of the ground state is now 
a scalar, and the atomic dipole operators b ± (r) [Eq.(19)] 
reduce to c-number functions that are simply given by 
the electric field profile 



b-(r)=£(r) 



(60) 



The advantage of such a transition is that it provides a 
good way to single out the effect of the interface on the 
basic atomic external dynamics, with the minimum com- 
plications introduced by the interface-modified internal 
atomic dynamics. 



1. Fluorescence rate 

The atomic ground state reducing to a single level, the 
only nontrivial feature of the optical pumping equation 
( |33| ) is the total atomic fluorescence rate r'(z), that is 
given by the trace of either term of the right-hand side 
of Eq.(M). One thus finds 



r>(z) = r' 00 (c ll (z)\z oll \ 2 + c x (z)\to±\ 2 )e 



(61) 



with £0 = (£o||j£o-l)- The influence of the interface on the 
internal ground state dynamics hence amounts to a dif- 
ferent space-dependence of the broadening of the ground 
state level in addition to the one associated with the 
space dependence of the driving laser field. For compar- 
ison, the parenthesis in Eq.(pl|) is space independent in 



free space and equal to |£o| 2 - It is also interesting to note 
that the fluorescence rate ( |6l| ) involves the intensity of 
the 'transverse' (|£o_l| 2 ) and 'longitudinal' (|£o||| 2 ) parts 
of the driving field independently, as a result from the 
rotational symmetry breaking due to the interface. As 
a consequence, the two elementary polarizations of the 
evanescent wave yield different spatial variations of the 
fluorescence rate: 



r'^)(z) = r' ooC ||(z)e- 



■2kz 



v(TJVf) 



k 2 Cii (z) + Q 2 c ± (z) ^_ 2kz 



(62a) 
(62b) 



Hence, depending on the polarization of the driving 
evanescent wave, the atomic fluorescence permits to 
probe different combinations of the correlation tensor 
components cu(z) and c\_(z). 



2. Radiation pressure force 



By inserting the field correlation tensor ( |C19 , C23| ) 
and the evanescent field profile (|5^) into the general re- 
sults ( B9l P01 ) , one readily finds that the departure contri- 
bution^!]]) to the radiation pressure vanishes for a scalar 
atom. After some algebra, the feeding term (|9|) yields 
the following result for the radiation pressure force 



F>P>(r) =T'(z)%Qe x + 

+2r' oc 1%k ai(z; 0) Im (£ ||£oj 



-2kz 



(63) 



The first term of the force ( |63| ) corresponds to the rule- 
of-the-thumb expression for the radiation pressure: it is 
the product of the fluorescence rate (61) and the photon 
momentum HQ carried by the evanescent wave along its 
propagation direction. The second term arises from the 
fact that the actual driving field consists of the sum of the 
incoming laser evanescent wave and of the reflected part 
of the field radiated by the atomic dipole. Because the 
radiation pressure force exerted on a J g = — » J e = 1 
atom is proportional to the phase-gradient of the total 
driving field, it actually appears as the sum of a contri- 
bution proportional to the phase gradient of the evanes- 
cent incoming wave [first term of Eq.(|63|)l and of a term 
involving the phase gradient of the reflected dipole field 
[second term of Eq.j63|)]. It is therefore not surprising to 
find that the reflected field contribution to the radiation 
pressure is proportional to the axial coefficient a± that, 
as previously discussed, is directly connected to field re- 
flection processes at the interface. 
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Fi (propagation vector) 




FIG. 3. Illustration of the radiation pressure force in a cir- 
cularly polarized evanescent wave. 

Fi: usual radiation pressure (parallel to driving field's propa- 
gation vector); F2: correction due to partial field reflection at 
the dielectric surface; F: full force. The helicity vector h of the 
circular polarization and the Poynting vector of the evanes- 
cent wave (using the standard definition) are also shown. 

For the TE and TM polarizations of the evanes- 
cent driving field, this correction to the radiation pres- 
sure force is difficult to observe: indeed, the vector 
Im (£o||£o_l) van ishes in the TE-case and is parallel to 
e x in the TM-case \cf. Eq.(|59|)], thus modifying slightly 
the magnitude of F^ sp \ A more prominent modification 
occurs for a generic combination of TE and TM polar- 
izations. The effect is actually maximum for a circular 
polarization of the evanescent wave in a plane perpendic- 
ular to the interface, which can be achieved for 

^ ) =^™ ) +^r ) . (64) 

The field's polarization is then a + with respect to an axis 
parallel to the interface (given by the 'helicity' vector h = 
(2Q/k)[e x — (K/k)e y ], cf. Fig. |^). The total fluorescence 
rate is given by 

r ' (CT) W=T^^ (c {] (z)+c ± (z)) e" 2 - (65) 

and the reflected field contribution to the radiation pres- 
sure force is found perpendicular to the helicity h, yield- 
ing 

F( s p'(z) = %QT' {rT \z)e x 

-2?igr^ a x (z; 0) e~ 2KZ (te x + e„) . (66) 

Thus, the radiation pressure force forms a nonzero, z - 
dependent angle with the evanescent wave propagation 
vector Q e x . This is represented in Fig. || where the angle 
and magnitude of the radiation pressure force ([36]) are 
plotted as a function of the distance z from the dielectric 
surface. 




k z 



1 2 3 4 5 6 
k z 

FIG. 4. Direction of the radiation pressure force in a circu- 
larly polarized evanescent wave, as a function of distance from 
the dielectric. The curve shows the angle between the actual 
force and the evanescent wave's propagation vector (the vec- 
tors F and Fi of Fig. ^). For these parameters, the field's 
Poynting vector forms an angle of —45° with the propagation 
vector Fi. The inset gives the magnitude of the force, using 
either the exact result (^) (full line) or ignoring both the in- 
creased fluorescence rate and the force correction F2 (dashed 
line). 

Parameters: refractive index no = 1.5, evanescent driving 
field with k = k, Q — \/~2k. The force is plotted is units of 
TikT'^, as a function of kz. 

We finally note that a related effect was discussed by 
Roosen and Imbert [Jl3): these authors pointed out that 
the Poynting vector of an evanescent wave is not parallel 
to its propagation vector if the wave is circularly polar- 
ized. If one assumes the radiation pressure force to be 
parallel to the Poynting vector of the driving field, which 
seems tempting because the Poynting vector represents 
the local momentum of the field, one thus finds a result 
reminding of (|66|). In fact, it is well-known that this as- 
sumption is not correct, as can be readily checked on the 
simple example of a J g — — > J e = 1 atom interacting 
in free space with two plane waves 

£(r) = e lkx e z + e^e* (67) 

leading to a Poynting vector 

II oc e x + e y — cos [k(x — y)] e z (68) 

whereas the radiation pressure force is clearly oriented 
along e x + e y . This problem can be readily solved by 
noticing that the Poynting vector II is only defined up 
to a curl, as shown by the continuity equation 

V.II + gV = (69) 



where e is the electromagnetic energy density. It is thus 
straightforward to show that it is always possible to de- 
fine a "novel" Poynting vector satisfying the continu- 
ity equation and being parallel to the radiation pressure 
force. 
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3. Momentum diffusion tensor 

The momentum diffusion tensor can be readily ob- 
tained from Eq.(45) using the field correlations (C1E, 
C23| ) and replacing the dipole operators b ± (r) by the 
evanescent field profile according to Eq.(|58|). One thus 
finds 



D 



8 oz z 



h 2 r< 



d 2 A 



2 ds k dsi 



[z;0) (70) 



where 

A(z; s) = & ' j (z; s) fa £ oj exp(-2«z + iQ s.e x ) (71) 

We now discuss the physical significance of Eq.(|7(]) in 
more details, with an emphasis on the influence of the 
interface on momentum diffusion. We start by consider- 
ing the first term on the right-hand side of Eq. (|70| ) , that 
corresponds to depart t see Eq.(^g)]. In order to single 
out the interface contribution, it is convenient to write 
the total fluorescence rate in the form 



2 £ -2 KZ 



+ T 'int{ z ) 



yielding two contributions to -Depart 



D 



depart 



# d 2 T' mt (z) 
8 dz 2 



(72) 



(73) 



The first term on the right-hand side of Eq.(|73|) does 
not involve any surface-induced effects (apart from the 
existence of the evanescent driving field) and is associ- 
ated with the shrinking of the atomic spatial coherence 
(hence a broadening in momentum space) due to the 
non-uniform, exponential photon absorption probability. 
This corresponds to the intuitive fact that the coherence 
of an atomic wavepacket incident on an evanescent wave 
mirror will be destroyed more efficiently by absorption- 
spontaneous emission cycles close to the interface than 
far away as a result of the inhomogeneous laser intensity. 
The second term of Eq.(|7|) is a correction to -Depart 
arising from the modification of the fluorescence rate by 
the interface, that is responsible for an additional spatial 
modulation of the photon absorption probability, hence 
an additional cause for momentum diffusion. As is well- 
known, the dipole damping rates T||(z) and Tj_(z) dis- 
play an exponential dependence close to the interface, 
and then tend toward their asymptotic value with 
some damped oscillations ^J. Y' int (z) 1 and consequently 
^depart' are therefore expected to exhibit the same kind 
of behaviour [see Eqs. ( |6~t| ) and ((73|)]. 

Consider now the second term on the right-hand side of 
Eq.@, corresponding to D^ d [see Eq.@]. Again, the 
influence of the interface can be sorted out by expressing 
A in the form 

A(z;s)=A 00 (z;s)+A int (z;s ll ) (74) 

Aoo(z;s) = C%>(s)fa£oj exp(-2Kz + iQs.e x ) (75) 

Ain t (z;s\\) = C£? t (2;S||)f5i£<y exp(-2Kz + iQ s.e x ) (76) 



Because Aoo only involves the free-space correlation ten- 
sor C^, its contribution to i? fe ' ed is analogous to that 
usually encountered in free space for any driving laser 
field, i.e., it accounts for the random walk of the atoms 
in momentum space due to their recoil during sponta- 
neous emission of photons, the spontaneous emission di- 
agram being assumed as in free space. The contribution 
of Ai n t to momentum diffusion, involving the complex 
nonlocal correlations between the vacuum field compo- 
nents induced by the interface, exhibits an interesting 
feature. Because Amt is independent of the z-component 
of the relative position s (because such is C\^ t ), the 
interface-induced modifications of momentum diffusion 
due to changes in spontaneous emission in the vicinity of 
the dielectric medium will only take place parallel to the 
interface, a result that is not obvious when considering 
the important modifications of the spontaneous emission 
diagrams due to the interface Gj . 
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FIG. 5. Momentum diffusion coefficient for a scalar atom 
driven by a TE [(a)] and a TM [(b)] polarized evanes- 
cent wave. Thick solid line: trace D(z) ([7?]) of the 
diffusion tensor; dashed line: lateral diffusion coefficient 
Dn{z) — D xx (z) + Dyy(z); dotted line: diffusion coefficient 
D zz (z) perpendicular to the interface; thin solid line: trace 
D(z) of the diffusion tensor in free space (both the modified 
vacuum correlations and the curvature of the driving field are 
neglected) . 

Parameters: same as Fig. ^. The diffusion coefficient is plot- 
ted in units of h 2 k 2 T', as a function of kz. 
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Quantitative results for the momentum diffusion tensor 
are displayed in Figs. ||-|| In Fig. [| is plotted the trace 
of the diffusion tensor 

D{z)=Y J D i ' i {*) (77) 

i 

that permits to estimate the total width Ap of the atomic 
momentum distribution: for a spatially constant diffu- 
sion coefficient, 

Ap 2 ~ 2Dt (78) 
which may be generalized to 

Ap 2 ~ 2 J dtD[(z(t))} (79) 

where (z(t)) is the mean atomic trajectory in the evanes- 
cent field's dipole potential. Eq.(f79j) is valid provided the 
momentum diffusion is sufficiently small so that individ- 
ual atomic trajectories remain close to the mean path, an 
assumption that may become questionable depending on 
the experimental conditions. The integral d79) may be 
estimated from the typical timescale r ~ np z/ i nc for 
the reflection HH ( ~Pz,inc- incident atomic momen- 
tum along z) and the position zq of the turning point 
of the mean path: Ap 2 ~ 2tD(zq). The diffusion co- 
efficients plotted in Fig. || may thus be viewed as the 
squared momentum width of the reflected atoms, given 
the interaction time and varying the distance z = zq of 
the turning point. 

From Fig. ^, we observe that the atomic momentum 
diffusion perpendicular to the surface (the coefficient 
D zz represented by the dotted curve) is, on its own, com- 
parable to the diffusion in free space (the thin solid curve, 
neglecting the curvature of the exponential fluorescence 
rate and the modified vacuum correlations). This may 
be compared with Fig. |^ where the evanescent driving 
field has a large decay length 1/k (total internal reflec- 
tion close to the critical angle). The exponential field 
profile is then essentially constant and the perpendicular 
momentum diffusion decreases. In this direction, the mo- 
mentum diffusion coefficient D z ' z is now determined, on 
the one hand, by the spontaneous photons' recoil, D^ d , 
and, on the other hand, by the derivatives of the fluores- 
cence rates T^j_(z) that enter into -Depart (©■ Also, for 
the long-range evanescent wave, the TE and TM polar- 
izations produce a similar momentum diffusion since they 
both correspond to a (nearly) linearly polarized driving 
field. 
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FIG. 6. Same as Fig. |f| but for an evanescent driving field 
with larger decay length: k = 0.1 k, Q w 1.005 k. 

In conclusion, we observe that the momentum diffusion 
tensor is systematically larger close to the dielectric sur- 
face compared to free space (the thin solid lines in Figs. 0, 
|), by a factor of about three to four. Translating the 
width of the atomic momentum distribution into a coher- 
ence length, we see that close to the surface, spontaneous 
emission destroys the atomic coherence more efficiently 
than expected from free-space considerations. This may 
have been expected from the enhanced fluorescence rates 
rii x( z ) compared to Too, as well as from the spatial sub- 
wavelength structure of both the evanescent field and the 
vacuum field correlations. Finally, we note that this diffu- 
sion tensor also translates into an increased temperature 
limit for radiative atom traps in the vicinity of surfaces 




C. J g = 1/2 atom 

In this section, we turn to a second application of 
the theory developed in Sec. || and focus on the optical 
pumping processes taking place in atoms having a non- 
trivial Zeeman sublevel structure in their ground state. 
By considering the simple case of a J g = 1/2 atom, we 
show that these processes are modified both quantita- 
tively and qualitatively in the vicinity of the dielectric 
interface: since the fluorescence rates depend on both 
the atom-interface distance and the atomic dipole ori- 
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entation, the optical pumping rates are modified. More 
strikingly, a net ground-state magnetization is predicted 
to occur in a linearly polarized driving field if the atomic 
recoil is taken into account. We illustrate this effect by 
an explicit calculation of the radiation pressure force. 



1. Atomic magnetization variables 

In the case of a J g = 1/2 atom, the Wigner function 
W(r, p) takes the form of a hermitian 2x2 matrix de- 
scribing the populations and Zeeman coherences of the 
ground state. It is convenient to represent this matrix 



using the vector <x = (a x , a y ,a z ) of Pauli matrices 



W 



(w 



(80) 



With this definition, the scalar function w(r, p) describes 
the phase-space distribution of the total population, 
whereas the real vector J(r, p) (analogous to the Bloch 
vector for a two-level atom) gives the phase-space dis- 
tribution of the atomic magnetization: for unpolarizcd 
atoms, e.g., one has J = 0; for an ensemble of atoms 
prepared in the sublevel | + 1/2) z with respect to the z 
axis, e.g., J z (r,p) = w(r, p) and J x . y = 0. The evolution 
equations for the total population and the magnetization 
vector are obtained from the G.O.B.E. by taking the ap- 
propriate traces: 



w(r,p) = Tr JT(r,p), 
J(r, p) = Tr [ctW{v, p) 



(81a) 
(81b) 



We next need the action of the reduced dipole oper- 
ators b ± (r) on the Wigner matrix. From the Clebsch- 



Gordan coefficients for the J g — 1/2 — * J e 
transitions, one finds 



Jg, Jg 



b - (r) = 0£(T)+iatr x £(r), 

J e = 1/2 : 13= 1/3, 
J e = 3/2 : p = 2/3, 



a 



-1/3, 
1/3; 



1 

(82) 
(83) 



while b + (r) is given by the hermitian conjugate. The first 
term (3 £(r) in Eq.([32|) is similar to the reduced dipole op- 
erator for a scalar atom (it is parallel to the driving field's 
polarization vector £(r)); the second term ia er x £(r) ac- 
counts for couplings between the ground-state Zeeman 
sublevels, and therefore describes specific multilevel ef- 
fects such as Raman couplings^] 



l We refer to the paper by Dubetsky and Berman ^1 f° r 
a generalization of the decomposition ( B^| o f the dipole op- 
erator to the case J g > 1. Cf. also Ref. 
therein. 



and references 



2. Internal dynamics 

As in the scalar atom situation, we start our analysi s by 
considering the classical optical pumping equation ( |33| ) 
accounting for the internal atomic dynamics. In the 
J g = 1/2 case, two quantities characterize the internal 



atomic dynamics: the ground state light shifts [Eq.(f24|)] 
and the optical pumping rates [Eq.(^)] (we do not con- 
sider in this section the energy level shifts induced by 
the interface). We discuss in some detail the case of a 
circularly polarized evanescent wave and calculate the 
pumping rate. 

a. General. Let us first write down the light-shift 
Hamiltonian (|24|) (also a hermitian 2x2 matrix): 

H eff = hA'e- 2KZ b+ •£„ 

= HA' e~ 2KZ ([3\£ \ 2 + ah • er) (84a) 
h = Im£*x£ (84b) 

where bg is given by Eq. (j32[) , replacing £(r) with £ , 
and where the h vector corresponds to the helicity of the 
driving field. Following Cohen- Tannoudji and Dupont- 
Roc pH , one may interpret the second term of the light- 
shift operator (p4|) in terms of a fictitious magnetic field 
parallel to the helicity h. This interpretation is supported 
by the equation of motion for the atomic magnetization 
vector J: it is obtained using Eq.(|34|) for the Wigner 
matrix W that takes the following form 



W r . 



r> oo e- 2KZ C i > i (z;0)(buWb+ 



(85) 



After some straightforward algebra with the Pauli ma- 
trices, we find from the light-shift Hamiltonian ( pi) ) and 
Eq. (p5[) the following equation of motion 



dJ_ 

~dt 



= 2a A'e- 2KZ h x J + 2a 2 r^ e~ 2KZ C(z)h 



({Ty T)C{z)J + {TrC{z)) T3 
+ 2T' oc e- 2KZ (a 2 {C{z), T}J - a(3[C(z), jF]j) 
We have used the following tensors 



{C{z)f 3 = C n 'Hz;0) 



C||(z) 
C||(z) 
c±{z) 



(86b) 



x,y.z 



W=Re&Gy (86c) 
Since the total population w = Tr W is conserved by 



optical pumping, we have put w = 1 in Eq. 

In the equation of motion (86a), one identifies the 
precession of the magnetization vector J around the ef- 
fective magnetic field vector h (the first term) and the 



14 



feeding of the magnetization J through optical pump- 
ing (the second term). The third and fourth terms de- 
scribe the damping of the atomic magnetization through 
absorption-spontaneous emission cycles. The effect of the 
dielectric interface is encoded in the (diagonal) tensor 
C(z) whose elements are proportional to the dipole damp- 
ing rates Tu and T± [see Eq.(|57j)]. It is interesting to note 
from Eq.(|86|) that optical pumping creates a magnetiza- 
tion aligned parallel to the vector C(z)h that is generally 
not parallel to the helicity h as a consequence of the 
anisotropic fluorescence rates. For the pumping process 
close to the dielectric interface, the one-point correlation 
tensor C(z) hence plays the role of an (anisotropic) 'effec- 
tive magnetic susceptibility', linking the induced atomic 
magnetization to the effective magnetic field. 

b. Discussion of element ary p olarizations. According 
to th e ligh t-shift operator (84a) and the equation of mo- 
tion (36a), the optical pumping process is characterized 
by the helicity vector h of the driving field. For the ele- 
mentary polarizations of the evanescent wave, it becomes 



(a) TE polarization 




(b) circular polarization 




(c) TM polarization 




TE polarization: 
TM polarization: 



h (TE) 
h (TM) 



-(2nQ/k 2 )e y 



(87a) 
(87b) 



We also mention for later use the case of a circularly 
polarize d evan escent wave. As discussed at the end of 
Subsec. Ill B 2 , the field's polarization vector is given by 
Eq.(|64|) and one has 

circular polarization: = (2Q/k)[e x — (n/k)e y ] 

(87c) 

In the TE case, the helicity vanishes and hence no net 
magnetization builds up. The light-shifts for the two 
ground state Zeeman sublevels being identical (for the 
same reason), one would expect this case to yield a sit- 
uation similar to the scalar atom of Subsec. IIIB. Wc 



shall see, however, that this is no longer true when the 
atomic recoil is taken into account (cf. Subsec. [II C 3). 

In the cases of TM and circular polarization, the helic- 
ity is nonzero and optical pumping leads to a net atomic 
magnetization. In steady state, it aligns parallel to the 
y-axis in the TM case, and parallel to an axis in the xy- 
plane for a circularly polarized evanescent wave. In both 
cases, the light-shift operator is diagonal with respect to 
these axis', and the anisotropic magnetic susceptibility 
does not come into play. [This would be the case for a 
different relative phase between the TE and TM polar- 
izations in Eq.(J64|), giving the helicity vector a nonzero 
component h z ] 

c. Example: circular polarization. As an example, 
we study in more detail the case of a circularly polarized 
evanescent wave and defer the TM case to Appendix [j^. 



FIG. 7. Illustration of the fluorescence cycles for a 
J = 1/2 — > J e —3/2 atom driven by (a) a linearly polar- 
ized field (TE polarization); (b) a circularly polarized field 
(combination of TE and TM polarization); (c) a TM polar- 
ized field (elliptic polarization). 

For clarity, only the fluorescence cycles starting from the sub- 
level | — 1/2} are shown. Note that these sublevels are defined 
with respect to a quantization axis that, in case (b), differs 
from cases (a) and (c). 

It is convenient to introduce a rotated coordinate sys- 
tem with the e x i unit vector being parallel to the helic- 
ity (37c). Choosing this as the quantization axis diag- 
onalizes in fact the light-shift operator (p^). From the 
magnetization's evolution Eq.(|8^), it may also be veri- 
fied that the components J y i , J z decouple from J x > . If 
we assume that the atoms are initially unpolarized, the 
pumping process only depends on J x i whose evolution is 
given by 



dJ x > 
dt 



= T p (z)[l- J x >] 



where the 'optical pumping rate' equals 

>.Q 2 



r P («) 



4a z — T' e 

£2 00 



-2kz 



C||0) 



(88a) 



(88b) 



The significance of this equation becomes evident if we 
write down the evolution of the populations w± = ^(1 ± 
J x i) of the Zeeman sublevels | ± 1/2) x < with respect to 
the quantization axis. The following rate equations are 
easily found 



dw+ 



r p (z)w_, 



= -T p (z)w- 



(89) 



and we see that the pumping rate T p (z) governs the tran- 
sition | — l/2) x r — > | + l/2) x ' between the Zeeman sub- 
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levels. In this transition, the atom absorbs a <j + polar- 
ized photon (with respect to the x' axis) from the driving 
field and spontaneously emits a 7r polarized photon (c/. 
Fig. 0b). The reverse transition is impossible since the 
driving field's polarization is purely a + . From this el- 
ementary pict ure, w e may understand why the optical 
pumping rate (88b) involves the coefficient cii from the 
vacuum correlation tensor: the spontaneous n photon 
has in fact an electric field parallel to the x' axis, and 
its emission rate is proportional to the strength of the 
vacuum fluctuations polarized along this axis, hence pro- 
portional to the element C x ' x = c\\ of the correlation 
tensor. 

The rate equations ( |S9| ) show that in steady-state, the 
atoms are completely magnetized along the helicity vec- 
tor of the driving field (the x' axis). With respect to op- 
tical pumping in free space, the only difference is hence 
the nonzero angle between this vector and the evanescent 
wave's propagation vector Qe x . 

A more detailed investigation of the pumping process 
may be done in the transient regime where the interaction 
time t is smaller than the pumping time 1/T p . This 
regime is in fact typical for atomic mirror experiments 
where one seeks to avoid spontaneous emission because it 
reduces the coherence of the reflection. An approximate 
solution of Eq.(jS8|) in the transient regime is (for initially 
unpolarizcd atoms) 



T p t < 1 



4 a 2 ^r^e- 2KZO c,|(z )r (90) 



J x > = w + -w 

>.Q 2 
k 



The population difference now depends on the coefficient 
cm (zq) at roughly the distance of closest approach zq. The 
estimate (^0|) is actually very crude, since it neglects the 
fact that the sublevels | ± 1/2) x > are subject to differ- 
ent light shift potentials in the circular polarization case. 
This leads to a different potential (and, ultimately, ki- 
netic) energy after the sublevel change — a feature that 
has been studied already for both spontaneous |jj| , [jl|[l6| 
and stimulated |55|-]57|] transitions between sublevels. In 
order to describe both the center-of-mass motion and the 
anisotropic vacuum correlations, one may use the full 
Fokker-Planck equation derived in Sec. 0. An exam- 
ple of the corresponding 'recoil-induced magnetizations' 
is given in the next subsection. 

For the sake of completeness, we have analyzed in Ap- 
pendix [d] the optical pumping in a TM-polarized evanes- 
cent wave. In this case, the atoms do not become com- 
pletely magnetized because the driving field [Eq.(|59|)] is 
elliptically polarized. [It contains both o~ + and a~ com- 
ponents with respect to the y-axis, cf. Fig. [j]c.] We find a 
steady-state magnetization Jy Stot ^ = —2kQ/(k 2 + Q 2 ) > 
— 1, and a pumping rate that is again proportional to the 
coefficient cm . 

Summarizing, if atomic recoil is neglected, the internal 
dynamics close to the dielectric is subject to the follow- 
ing modifications as compared to free space. The optical 



pumping rates increase and differ according to the polar- 
ization of the spontaneous photon emitted in the pump- 
ing cycle. In a circularly polarized, evanescent driving 
field, a net atomic magnetization builds up that does 
not align parallel to the field propagation vector for two 
reasons: first, the magnetization is determined by the 
field helicity that is not parallel to its wave vector, and 
second, atomic magnetization and field helicity are con- 
nected by an anisotropic effective susceptibility because 
of the anisotropic fluorescence rates. 



3. External dynamics: recoil-induced magnetization 

We now consider the G.O.B.E. accounting for the 
atomic recoils during absorption and emission of pho- 
tons. Among the huge diversity of non standard effects 
expected in the external dynamics of multilevel atoms in 
the vicinity of a vacuum-dielectric interface, that could 
not be addressed in a single paper, we focus here on ra- 
diation pressure and show that this force may induce a 
net magnetization of the atomic ground state for certain 
classes of the atomic velocity in a situation where the 
classical optical Bloch equations would predict a zero re- 
sult (linearly polarized driving field). 

a. Radiation pressure force. Consider the situation 
of a J g = 1/2 atom driven by a TE polarized evanes- 



cent wave. Using the general expressions ( p9| , 41) for 
the radiation pressure force, after a straightforward cal- 
culation using the Pauli matrices, one finds that the de- 
parture contribution F^p art ( |4l"| ) vanishes, whereas the 

feeding contribution (39) is a sum of two terms 

analogous to the ones encountered for the scalar atom 
[Eq.([33|)1. The first of these contributions, , corre- 
sponds to the product of the fluorescence rate and the 
evanescent field phase gradient. This force is hence par- 
allel to the evanescent wave propagation vector Qe x . We 
find that this term does not couple the population w 
and the magnetization components J. The second con- 
tribution to the radiation pressure operator (|39|), , 
is associated with the gradient of the field emitted by the 
atomic dipole and backreflected towards the atom by the 
dielectric interface, and involves the axial part of the field 
correlation tensor. This contribution, that vanished for 
a TE polarization and a J g = — * J e = 1 atom, takes 
a nonzero value in the present situation and gives rise to 
a magnetization-population coupling. We interpret this 
coupling in terms of generalized rate equations for the 
Zeeman sublevel phase-space distributions. 

To be more explicit, the contribution F^ to the radi- 
ation pressure force appears in the following way in the 
Liouville equatio n for the magnetization component Ji 
[cf. Eqs.(p5[[81b|)l: 



Tr {a, (r) • V p W} = F\$> (r; J,) • V p J, 



( S P); 



(91) 



1G 



It is characteristic for the force F;*?^ that on the rhs 
only Ji appears. A similar result holds for the total pop- 
ulation w. The following expressions are found (the z- 
dependence of cm,c_l has been suppressed for clarity): 

Fj^(r;u;) = F^' 1 ) [/? 2 C || + a 2 ( C ,| + c±)] (92a) 

FjJ (r; J x ) = F^ 1 ) [/3 2 C || - a 2 ^ c ± )} (92b) 

F<^(r; J y ) = F^' 1 ) [/3 2 Cjl - a 2 (c,| + Cj _)] (92c) 

Fj^(r; J,) = F^' 1 ) [/3 2 C |, + a 2 (cj| - Cj _)] (92d) 



F (ap,l) = p/ e -2«^Q ea 



(92e) 



Note that the effect of the radiation pressure force ( |92j ) 
differs between the magnetization components. This fea- 
ture is interpreted below using the rate equations j95| ) 
where we show that the radiation pressure depends on 
whether or not a fluorescence cycle leads to a sublevel 
change. In any case, the contribution ( p2|) to the radi- 
ation pressure is parallel to the propagation vector Qe x 
of the evanescent wave. We also observe that the effect 



of the force (92a) on the population w is proportional to 
the total fluorescence rate r'( TB '(z): 



p(sp) 

(1) 



(r;w)=hQe x r'( TE \z) 

= TiQe x (L'(z) + r'(z)) 
T' 7r (z)=(3 2 T' O0 e-^c ll (z) 
T' (7 (z)^2a 2 T' 00 e^\{c,(z) + 



(93a) 
(93b) 
(93c) 



The notations T'^ and Y' a refer to a quantization axis 
chosen along the y axis (cf. Fig. ^a): with respect to 
this axis, the (linearly polarized) driving field excites a 7r 
transition and T' v gives the fluorescence rate for sponta- 
neous photons with tt polarization (electric field parallel 
to the y axis). As discussed in the example of a circular 
driving field, this rate is proportional to the coefficient 
en . We observe that the fluorescence rate for a polar- 
ized photons is proportional to |(cii + cj_), these photons 
having an electric field in the xz plane. 

The contribution F^? to the radiation pressure oper- 
ator (39), in contrast to Eq.(^), mixes the population w 
and the magnetization J. More precisely, the Liouville 
equations for w, J x and J y contain the following terms 
proportional to the axial coefficient ai of the field corre- 
lation tensor: 



Tr{Fgf(r) 



TrK(F^(r).V p W0} 
V P W)} 



V P W} = f 



- f(sp,2) 



a/3 



dJ x 



(2) 



d P y 

a(3f {sp ' 2 



O dJy 
t - 

dp x 



Tr{a y (Fg)(r) 



a 2 / (sp ' 2; 
2KZ hk ai (z;0) 



dw 
dp v 
dw 
dp x 



(94a) 

(94b) 

(94c) 
(94d) 



The J z magnetization component is not coupled to 
{w, J Xl J y } in this case. 



b. Rate equations. In order to make the physical con- 
tent of the Liouville equations (^2[ ^J) more transpar- 
ent, we again consider initially unpolarized atoms and 
suppose that their momentum distribution is uniform in 
the y direction (perpendicular to the evanescent wave's 
propagation vector). In these conditions, it is possible to 
neglect terms involving the derivatives d/dp y in Eq.(p4[). 
The coupled Liouville equations then transform into a 
pair of rate equations involving only the sublevel popu- 
lations w± = h(w ± J y ) with respect to the y-axis: 

( ^lo+i + |-'Vr)w±+ 7±^t( z ) w ± - 7t-±( z ) w t 



-F±_±(aO • V p w± + F T ^±(z) • V p w T = 



(95a) 



The different quantities in these equatio ns a re easily 
found by comparison between the Bloch (36a) and the 
Liouville equations (§2, w&) : 



7±-t( z ) 
F±^±(z) 

F=F-±(*) 



2/3 A' e - 2KZ h K e z + T'^(z)hQe x 
K(z)hQe x T a 2 f^e x 



(95b) 
(95c) 
(95d) 



The significance of these results is clear. The J±^^(z) 
are the transition rates for a sublevel change | ± 1/2) 9 —> 
|=Fl/2) y ; both transitions take place at the rate T' a (z), the 
fluorescence rate for a ± polarized spontaneous photons 
[cf. Eq.© and Fig.Ra]. 

The forces F±_»±(z] lare the sum of the dipole force 
(the first term in Eq.( |95cD ) and the radiation pressure 
force due to fluorescence cycles where the atoms fall back 
to the same initial sublevel (the second term) . Since the 
driving field is linearly polarized, this force is propor- 
tional to the fluorescence rate T' n (z) for tt polarized pho- 
tons (cf. Fig. @a). 

Finally, the forces F T ^± are radiation pressure forces 
due to sublevel-changing fluorescence cycles | =p l/2)„ — > 
| ± l/2) y . They differ in two respects from the previous 
force. First, their mean value (averaged over the sub- 
levels) is proportional to the emission rate F' a (z) for o± 
polar ized photons. Second and more striking, the forces 
(|95cj ) are not the same for the transitions | + 1/2) y — > 
| — 1/2) y and | — \/2) y — > | + l/2) y , their difference being 
proportional to the weight function ax(z; 0) for the axial 
part of the field correlation tensor [cf. Eq.( |94cj )]. To un- 
derstand this result, we recall that for a scalar atom, the 
axial part comes into play when the atom, driven by a 
circularly polarized field, emits circularly polarized pho- 
tons. More precisely, the photons must be polarized in a 
plane perpendicular to the interface (the y'z plane in the 
example studied in the preceding paragraph). If this is 
the case, the axial correlation tensor results in a force in 
the polarization plane and parallel to the interface, with 
a sign depending on the helicity of the spontaneous pho- 
ton (see Fig. ||). We encounter here a similar effect for a 
J g = 1/2 atom: even in a linearly polarized driving field, 
the spontaneous photon's polarization is indeed circular 
as soon as the atom changes sublevel. Consider for exam- 
ple the transition | — 1/2) y — > | + l/2) y shown in Fig. 0a. 
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The spontaneous photon is a~ polarized and since the 
quantization axis is the y axis, its electric field lies in the 
xz plane. The emission of this photon hence gives rise to 
a force correction parallel to e x . On the other hand, the 
reverse transition | + 1/2) y — * | — 1/2) y is associated with 
a a + polarized photon and a force correction of opposite 
sign. 




FIG. 8. Illustration of 'recoil-induced magnetization' in a 
TE polarized evanescent wave. The momentum distributions 
for the sublevels | ± 1 /2) y (denoted by '+' and '— ') are shifted 
by the radiation pressure forces F±. Due to the difference in 
recoil momentum per pumping cycle, the sublevel distribu- 
tions separate and merge periodically at the optical pumping 
rate. The dashed line shows the net magnetization J y (p x ) 
as a function of the momentum component p x parallel to the 
evanescent wave's propagation vector. 

c. Experimental signature. As a consequence of 
the difference between the radiation pressure forces 
F T ^±(z), the atomic Zeeman sublevels absorb different 
momenta per optical pumping time. Their distributions 
hence separate in momentum space. But since the sub- 
levels have been exchanged in the pumping cycle, the sub- 
level momentum distributions w±(p x ) merge again after 
a second pumping cycle. The process is then repeated 
periodically in time, as shown schematically in Fig. |^. 
The maximum sublevel separation in momentum space 
is of the order of a few photon momenta 



max Sp 



F' 



-4?ifce, 



ai(z;0) 



' c\\(z) + cj_(z) 



(96) 



We note that in the transient regime, this phenomenon 
may be observed experimentally using similar techniques 
as for atomic diffraction experiments at normal incidence 
p8[ . For longer interaction times, one could think of Ra- 
man spectroscopy to detect the Zeeman sub leve l imbal- 
ance as a function of the velocity p x /M: Eq.(Mc) indeed 
predicts a magnetization u>+ — w_ = J y proportional to 
the derivative dw/dp x of the atomic momentum distri- 
bution ('recoil- induced magnetization'). Note that the 
periodic separation of the sublevels might be difficult to 
observe because their momentum distributions are also 
broadened by spontaneous emission (momentum diffu- 
sion) . 



IV. CONCLUSION 

We have formulated a theoretical framework to de- 
scribe the motion of an atom that is fluorescing in an 
environment with modified electromagnetic field modes. 
We provided expressions for the radiation pressure force 
and the momentum diffusion tensor in the limits of low, 
but semiclassical velocites and low saturation. This gen- 
eral theory applies to atoms with arbitrary Zeeman sub- 
level structure and environments with arbitrary elec- 
tromagnetic field correlations. One important result is 
that the internal and external (center-of-mass) dynam- 
ics of the atoms is determined by the two-point correla- 
tion tensor of the vacuum field around the atomic posi- 
tion. We then made explicit predictions for simple atoms 
(J g = 0, 1/2) in the vicinity of a flat dielectric surface. 
These results allow for precise estimates of the effects 
of spontaneous emission when atoms are either reflected 
from an evanescent wave mirror [ p^ , pJ3||lS| ] or trapped in a 
two-dimensional waveguide-like field configuration in the 
vicinity of a surface J19|-p5[. 

Even for a scalar atom (J g = ground state), the 
radiation pressure force exhibits quantitative and quali- 
tative changes with respect to free space. It is increased 
both due to the subwavelength structure of the evanes- 
cent driving field and the modified vacuum correlations. 
In particular, the radiation pressure is no longer paral- 
lel to the phase gradient of the driving field due to the 
partial reflection at the dielectric interface. The optical 
pumping of a J g — 1/2 atom in an evanescent wave shows 
similar modifications. The atomic magnetization vector 
is related to the field's helicity in an anisotropic manner, 
and even a linearly polarized field gives rise to an imbal- 
ance of the atomic sublevel populations in velocity space 
('recoil-induced magnetization'). The sublevel-selective 
detection of atoms reflected from an evanescent wave mir- 
ror is thus a sensitive probe of the electromagnetic field in 
the simple half-cavity realized by the vacuum-dielectric 
interface. 

The present work might be pursued in two directions: 
first, one could further explore the properties of radiation 
pressure in evanescent waves and consider situations be- 
yond the simple models studied in this paper. Let us 
mention some topics of particular interest: angular mo- 
menta J g > 1 because such atoms are actually used in 
experiments; the coupling between Zeeman sublevel pop- 
ulations and coherences for more complex field polariza- 
tions (TM, circular); and the momentum diffusion for 
Zeeman-degenerate atoms ]59[ | . One may also expect that 
the combination of relaxation processes and the center- 
of-mass motion in different light-shift potentials leads to 
a variety of motion-induced magnetizations, similar to 
the case of conservative couplings explored in reflection 
beam-splitters |35| and diffraction gratings |5f||57||. A 
second direction opens up if one considers different ge- 
ometries: cold atoms trapped in high-quality cavities, 
e.g., are currently receiving much interest in the fields 
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of cavity QED and nonlinear quantum optics. The gen- 
eralized optical Bloch equations derived in Sec. || could 
be used in their present form to study atomic motion in 
semiclassical driving fields and might be generalized to 
describe more complex phenomena as, e.g., the absorp- 
tion and transmission of a probe field or the influence of 
the atoms on the cavity properties. 
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APPENDIX A: DERIVATION OF THE G.O.B.E. 



i+Ai 



(15) 



We outline here the derivation of the quantum- 
mechanical master equation (|l5|). The only difference 
to the usual treatments J2{| is the quantization of the 
atomic center-of-mass motion, i.e., we take care of the 
ordering of the atomic position and momentum operators 
R, P. We only present the case without an external driv- 
ing field, since this field may be easily accounted for by 
adding a commutator with the interaction Hamiltonian 
Val (H) to the master equation of the reduced density 
matrix. We also assume that the electromagnetic field is 
at zero temperature (in the vacuum state) , as is usual for 
optical frequencies. 

In the interaction representation (with respect to the 
free atomic Hamiltonian Hq plus the free vacuum field 
Hamiltonian Hr), the evolution of the full atom + reser- 
voir density matrix par is given by 



PAR = [VAR(t),pAR] 



(Al) 



with the atom-field interaction given by the electric 
dipole interaction 



V AR (t) = -Vd(t)--E{R,t) 



(A2) 



Here, D(i) = Vd(t) is the atomic dipole operator and 
E(r, t) the electric field operator in the Heisenberg pic- 
ture (these operators evolve in time according to the 
Hamiltonian Hq + Hr). Anticipating the appr oximation 
of slowly moving atoms, we neglect in Eq.(A2) the time- 
dependence of the position operator R due to Ho (free 
flight). 



To solve Eq.(Al) in second-order perturbation theory, 



PA R(t + At) = p AR (t) + — I dr [V A r(t), Par(t)] 



(A3) 

This equation is iterated by inserting a similar expression 
f° r Par(t) under the integral sign. Taking then the trace 
over the vacuum field variables, one obtains the master 
equation for the reduced atomic density matrix p. As 
usual in perturbation theory, the second-order term in 
Var is simplified by factorizing the full density matrix 
according to par(t) — p(t) <8> p vac . The resulting master 
equation reads in the position representation 

t+At t' 

(r 1 |p(t + At)|r a ) = (r 1 |p(t)|r a ) + ^ j dr j dr E 
{ - (sW(r 1 y)£7 j ( - ) (n > r))d i (r')d j (r)(r 1 | /J (r)|r 2 ) 

-{E { +\v 2 ,T)Et-\v 2 y)){v l \ P {T)\V2)d l {T)d ] {T>) 

+ (E^(r 2 ,r)E ( -\r 1 ,T'))d 3 (r')(r 1 \p(r)\r2)d l (r) 

+ < J B| +) (r 2 ,r')^- ) (r 1 ,r))d J (r)<r 1 |p(r)|r 2 )rf l (r')}. (A4) 

For brevity, we did not write out the vacuum state for 
the field expectation value (0|^ +) (r x , T r )E { r } (r 2 , r)|0) 
where E^ (r, t) are the positive and negative frequency 
components of the field operator. 

We now explicit the free evolution of the atomic dipole 
operator for a two-manifold system 



d(r) = d- e - luJA ( T -^ +d+e lu ' A ^- t \ 



(A5) 



where d^ are the dipole raising and lowering operators 
in the Schrodinger picture at time t. We observe that 
the correlation time of the vacuum field fluctuations is 
much shorter than the timescale At for the evolution of 
the atomic density matrix. This imp lies that we may 
compute the time integrals in Eq.(A4) in the usual way 
p7[ : replace r by t in the argument of the density ma- 
trix and take the latter outside the integral; change to 
the integration variable r' — r and replace its border At 
by infinity; discard terms oscillating at twice the optical 
frequency; identify the Fourier transform of the two-time 
vacuum field correlations at the atomic frequency: 



£* J (ri,r 2 ) = 



(A6) 



= / d(T , -T)(E$ +) (r 1 ,T , )E! j -\r 2 ,T))e i »A ^' ~^ 



we first re-write it as an integral equation: 



The remaining time integral then turns out to be pro- 
portional to At, and after some term rearrangements, 
one obtains the following form for the master equation 
(summation over repeated indices is understood) 
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p(t + At)-p(t) 



At 



(A7) 



"El. 

2h 2 

V 2 



{v 1 \{£^{^B)dldT,p{t))\r 2 ) 



+ F ^(r 2l r 1 )dJ(r 1 |p(t)|r 2 )dt 
+ l<r 1 |[(/r A (R)- J ff ),p(t)]|r 2 ) 

We may identify the lhs of this equation with the time 
derivative (ri|p|r2) since the reduced density matrix 
evolves slowly on the time scale of the vacuum fluctu- 
ations. 



We finally recall that the correlation function (A6) is 
identical, up to a normalization, to the field correlation 
tensor defined in Eq.(p7|): (V 2 /h 2 )£^ = T^C^ . The 



first two lines of Eq.(A7) are then readily identified with 
the relaxation part p re iax of the G.O.B.E. (|l5|). Further- 



more, the last line of Eq.(A7) contains the level shifts 
(Lamb-shifts) due to the coupling to the reservoir. These 
shifts contain both the renormalization of the atomic 
Hamiltonian, Ha,oo — Hq, and its the interface-dependent 
part AH a (R) appearing in Eq. (||) . For simplicity, we do 
not write down explicit expressions and refer to Ref. ]4l| ] 
for a discussion of the atomic level shifts in the vicinity 
of a vacuum-dielectric interface. 



APPENDIX B: ADIABATIC ELIMINATION OF 
THE OPTICAL COHERENCES AND THE 
EXCITED STATE POPULATION 

In this appendix, we analyze in detail the validity con- 
ditions for the adiabatic elimination of the optical coher- 
ences and the excited state density matrix. 



1. Optical coherences 

The atoms are driven by a laser field that we describe 
by a monochromatic classical field 



£ L (r,t) = £ L {r)e- luJLt + c.c 



(Bl) 



The interaction Hamiltonian (||) becomes, in the rotating 
wave approximation, 

V AL = -V (d+.£ L (r)e- luJLt + d~ £* L (r)e iu,Lt ' 

= -V£ (d+ 4{r)e- luJLt + d~.£* (r)e Mit ) (B2) 

where the dimensionless vector £(r) for the field profile 
[Eq.(20)] has been used. The time-dependence of the 
interaction Hamiltonian is removed by passing into the 
"rotating frame", i.e., we write the optical coherence in 
the form 



Peg ^ Peg 



(B3) 



The equation of motion for p eg is now readily obtained 
from the Bloch equations 0, [§[ [l5]) and reads 



(B4) 



1 

ih 

1 

ih 



P e AH A (R)P e p eg - p eg P g AH A (R)P g 



2M 



Peg 



W£ Q 



([d+.£(R)]<r-p ee [d+.£(R)] 



The ordering of the terms takes into account that R is 
the atomic position operator. We have also introduced 
the abbreviation 



£ e (R) = C lJ (R,R)d+cn 



(B5) 



Eq.(B4) may be approximately solved if the detuning A 
is outweighing all the other frequencies. We therefore 
assume that the atoms are driven off-resonantly and at 
low saturation, |A| ^> r^, T>£ /fi, \AHa/U\, as in con- 
dition (fio"|). The frequency associated with the kinetic 
energy operator may be estimated in the Wigner repre- 
sentation [cf. Eq.(|35|)]. If we assume that the atomic po- 
sition distribution varies at most on the scale of the opti- 
cal wavelength, this term is overestimated by the Dopplcr 
shift kp/M. In this way, we find the third condition ap- 
pearing in @, |A| > kp/M. 

Given these conditions, the adiabatic solution to 



Eq.(B4), correct to first order in r M /A, reads 



D£q f Xoo /tjn 



[d+.£(R)]<7-p ee [d+£(R)] 



(B6) 



Note that the optical coherence p ge is equal to the her- 
mitian conjugate of (B6). We shall see in the next 
paragraph that the excited state density matrix p ee is 
muc h sm aller than the ground state density matrix a 
[Eq. (|B12[ )]. We may therefore neglect the former in 
Eq. (p36|) . Using the definition ( |l9|) of the reduced dipole 
operator b~ (R) and recalling that the projection opera- 
tor dfd~ — P e acts as the identity onto the excited 
state manifold, one sees that Eq.(B6) yields the expres- 
sion (p7|) for the optical coherences. 



2. Excited state 

The generalized optical Bloch equation for the excited 
state density matrix p ee reads 



i 

ih 



2M 



£ e (R), p e 



P e A£T A (R)P e , Pe 



W£ 



■d+4(R)]Pse ~ Peg[d-4* (R) 



(B7) 
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Inserting the adiabatic solution ( JB6| ) for the optical co- 
herences, one obtains 



Pee + £e(R), Pe 



2 

1 

it) 

Too so 
"~2~T 
Too s / , 



1 

if) 



2M 



, Pe 



SO 



P e AflA(R)P e - ftA-^[d + .£(R)][d-.^'(R)], Pe 



&(R), [d+ £(R)]cr[d-.€*(R)] 



v [d+.|(R)][d-.r(R)]Pee^e(R) 

+a e (R)Pee[d+.€(R)][d-.C*(R) 



(B8) 



We have used the saturation parameter s$ defined in 
Eq.©. 

To solve this equation approximately, we observe that 
the last term on the rhs is small compared to the term 
involving r m on the lhs, because of the low saturation 
limit so -C 1 [condition (|jo|)] . We also want to neglect the 
first two lines on the rhs (the kinetic and potential energy 
operators in the commutator). We noted above that the 
kinetic energy corresponds to a rate smaller than roughly 
the Doppler shift. It is hence negligible if condition (48), 
Too 3> kp/M, holds. We note that this condition may be 
re-written as 



(B9) 



i.e., the atoms move much less than a wavelength during 
the lifetime of the excited state. To estimate the potential 
energy term, we again use the Wigner representation and 
find [cf. Eq.(p5|)] that it is of the order of the force F e 
acting on the excited state divided by the width Ap of 
the atomic momentum distribution. We therefore need 
to suppose Too > F e /Ap [cf. Eq.(pl)] or 



< Ap 



(BIO) 



The momentum the atoms gain during the excited state's 
lifetime is hence negligible compared to the width of the 
momentum distribution. 

Given these conditions, we are left with the equation 



<?e(R)i Pe 



|[d + £(R)]a[d-.C(R)] 



(Bll) 



Using the fact that the spontaneous emission rates are 
positive for any polarization of the spo ntane ous photon, 
it is easy to prove that the solution to (Bll ) is given by 



Pe 



SO 



[d+.$(R)]<7[d-.r(R) 



(B12) 



3. Ground state 

The Bloch equation for the ground state density matrix 
a reads in the position representation 



(ri|d|r 2 ) = 4( r i 
in 



P 2 

2M 



P g AH A (K)P g , a 



r 2 



[d-.r(ri)]<ri|Pe fl |r 2 > - (r 1 |p 9e |r 2 >[d+.£(r 2 )]) 



-r oo C i J(r ai ri)d7(ri|/j ee |r2) 



(B13) 



We now insert the adiabatic expressions (BE , B12 ) for the 
optical coherences p eg and t he ex cited state density ma- 
trix p ee . The last line of Eq.( B13 ) readily yields the last 
line of the optical pumping equation ( p3| ) . In the second 
line, involving the optical coherences, we neglect terms 
of order Sq and obtain the following two contributions 

l(rx|[^ s (R), a] |r 2 ) - ^( ri \{g(R), a}\r 2 ) (B14) 

The commu tato r is due to the real part of the optical 
coherences (B6) (in phase with the driving field) and is 
characterized by the light-shift Hamiltonian 



ff,.(R) = ftA^[d-.e*(R)][d+e(R)] 
= ftA'[d-.r(R)][d+€(R)]. 



(B15) 



This Ham iltonian adds to the level shift P g AH A CR.)P g in 
Eq.( B13 ) to give the effective ground-state Ham iltonian 
H e ff(R) (p4|). The anticommutator in Eq.(B14) is due 
to the imaginary part of the optical coherences p eg (phase 
lag of order Too /A). It involves the ground-state operator 
G(R) defined in Eq.©: 



r' r 
-fg(K) = -fc* 

Too Sq 



2 2 



■*(R,R)&r(R)&7(R) 

[d-.^*(R)]g e (R)[d+.|(R)] (B16) 



We have thus obtained the optical pumping equa- 
tion (|23 



APPENDIX C: FIELD CORRELATIONS FOR 
THE VACUUM-DIELECTRIC INTERFACE 

In this appendix, we outline the calculation of the 
electromagnetic field correlation tensor for the vacuum- 
dielectric interface, following Carnaglia and Mandel fl48| . 
We assume that the dielectric fills the half-space z < 
and is characterized by the (real) refractive index tlq. 



This expression shows that the excited state density ma- 
trix is much smaller than that of the ground state, by a 
factor of the order of the sat urati on parameter sq. In the 
position representation, Eq.(B12) yields the result (28). 



1. Field modes 

Carnaglia and Mandel distinguish two types of electro- 
magnetic field modes for this geometry: 
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(a) modes incident from inside the dielectric and be- 
ing partially or totally reflected at the dielectric-vacuum 
interface. In the vacuum half-space, these modes are ei- 
ther propagating or evanescent, depending on the inter- 
nal angle of incidence. The wavevector of the incident 
wave in the dielectric is denoted ko-f, and the vacuum 
wavevector of the transmitted wave kf . Obviously, one 
has |koj| = no(uj/c) = ngk and |kf| = k. The wavevec- 
tors are decomposed according to kj = (kii,fe z ) where 
ky denotes the components parallel to the interface plane 
(the xy-plane) and k z the perpendicular component. The 
parallel components of ko and k coincide: k || = kii . The 
perpendicular components are such that in the dielectric, 
/coz > 0; in vacuum, k z is chosen such that k z > for 
the propagating modes and Im k z > for the evanescent 
waves. For later convenience, we introduce the abbre- 
viation u = |k|||/A;. Modes propagating (evanescent) in 
vacuum then correspond to < « < 1 (1 < u < no), 
respectively. 

We write f T (k 0T ,u;r) (fi = TE, TM) for the corre- 
sponding mode function that is normalized to unit (in- 
cident) amplitude in the dielectric. In the vacuum half- 
space, these modes have an amplitude equal to the Fres- 
nel transmission coefficient denoted by t(u, /i). Explicitly, 
one has (z > 0): 



f T (k ,M;r) - 
with k 



e x (u,ip,TE) 
e T («,^,TM) 

t{u,TE) 
t(u,TM) 



e^(u,ip,^)t(u,fj,) exp (ik T • r) , 
k(u cos ip, u simp, v), 
{0<u< n ), 

\/l-u 2 , 

(— simp, cost/?, 0), 

(v cos ip, v sin ip, —u), 

2^F^ 



(CI) 
(C2) 

(C3) 
(C4) 
(C5) 

(C6) 
(C7) 



Note that the TM polarization vector (C5) is complex 
for evanescent modes (u > 1). 

(b) modes propagating downwards from the upper 
half-space into the dielectric with wavevectors kj = 
(ky,— k z ) in vacuum and kox = (kn,— ko z ) in the di- 
electric. In the vacuum half-space, these modes con- 
tain a part reflected from the interface with wavevector 
kf. We note fj.(k^,//;r) the corresponding mode 



.0) 



function (normalized to unit incident amplitude in vac- 
uum) and r(u, fi) the Fresnel amplitude reflection coeffi- 
cient (z > 0): 



f J. (k, M5 r ) = exp (iky • ry ) (ej. (u, <p, fi) e 



- ik z z 



Jr) 



(u,ip,n)r(u,n)e 
with k| = k (u cos ip, u simp, — v) 
(0 < u < 1), 
ex (u, tp, TE) = (- sin ip, cos ip, 0) 



ik z z \ 



+ 

(C8) 
(C9) 

(C10) 



= e ( [\u,p,TE)^e^u,ip,TE) 
(u, cp, TM) = (v cos ip, v sin ip, u) (Cll) 
e[ (u, ip, TM) — (—v cos ip, —v sin ip, u) (C12) 
= -e 7 (u,ip,TM) 

V — \JtIq — u 2 



r{u,TE) 



r{u,TM) 



v + y/riQ — u 2 

n V - V n Q - 



(C13) 



(C14) 



UqV + y/nZ — u 2 

For a detailed discussion of the orthonormalization of 
these field modes, cf. Ref. (iq . 

2. Vacuum field correlation function 

Upon quantization, the electric field operator in the 
vacuum half-space may be written as a sum over the two 
types of modes introduced above, the mode functions be- 
ing multiplied by the usual creation and annihilation op- 
erators 48 6(J. Using the bosonic commutation rules, 
Carnaglia and Mandel obtain the following result for the 
vacuum correlation tensor (z\, z 2 > 0): 



|£! +) ( ri ,T)£H(r 2 ,0)|0) 



(C15) 



+ 



d 3 k ^ hu> 
(2^)3 2n 2 e 



d 3 ki Tilo 



/Ti( k oT, m; r i)/fj( k oT: m; r 2) z 



< (2tt)3 2e ^ 



/ii(kj.,M;ri)/|j(kj.,^;r 2 )e" 



The signs ">, <" on the integral signs are to remind 
that the wavevectors k f and kj only run through a half- 
space. 

From expression ( p!5 ), we may readily read off the 
Fourier transform £ * J (ri, r 2 ) of the correlation tensor 
defined in Eq.(A6), by using as integration variables the 
frequency u> and the polar coordinates u, ip of the in-plane 
vector kii/fc. After some algebra, one obtains the follow- 
ing representation of the normalized correlation tensor 



C^(n,r 2 ) = (C16) 
J dl P^ /ti( k 0T) M;ri)/fi( k 0T, M; r 2) + 



n 

3 f duu 



8tt 
3 

8^ 



,2 


2tt 



[From here on, the wavevectors k^j have magnitude 
|k+ 1 1 = loq/c — k where ljq is the atomic transition fre- 
quency] In the following, we show that the correlation 
function flClC| ) may be written as a sum of two parts, 
G = Coo + Ci n t , one corresponding to the free-space cor- 
relation function, and the other one representing the in- 
fluence of the interface. 
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a. Free-space part 

The free-space correlation function is obtained from 
those upward propagating modes that are homogeneous 
plane waves above the dielectric, on the one hand, and 
from either the incident or the reflected parts of the 
downward propagating modes, on the other. These con- 
tributions may be combined using the following property 
of the Fresnel coefficients 



1 



T t 2 {u,fi) +r 2 (u,n) = 1 (C17) 



u < 1 



that follows from the relation r = —r' where r' is the 
reflection coefficient for upward propagating modes (reci- 
procity), and energy conservation. One finally gets the 
result 



C^(rx,r 2 ) = |- JduJ 



dtp 



vT 



E 



(e Tl (u,^, M )e^-(u,(^)e lMzi Z2) + 

+ eii {u, <p, n)e* u (u, ip, n)e- lk ^- z ^ x 
x exp[tk|, • (rn ! - rn 2 )] (C18) 



If the integral is written in terms of the unit vector n = 
k/fc, one recovers the familiar expression for the free- 
space vacuum correlation tensor that appears in Eq.(18) 



Due to translational invariance of the vacuum field, the 
free-space correlation tensor only depends on the differ- 
ence vector s = r 2 — iv From rotational invariance, it 
follows that the tensor may be decomposed into a scalar 
part, proportional to the unit tensor, and a quadrupolar 
part with zero trace, proportional to s l s^ — ^s 2 6 h: * . More 
explicitly, the correlation tensor may be written for two 
neighboring points 

^'( S )=c^; ( S )+c^ 2 ( S ) 

w (l - ^fc 2 s 2 )<f <■> - ^fc 2 (sV - §s 2 S«) (C19) 

We have limited ourselves to second order in ks, which 
is sufficient to compute the radiation pressure force and 
the momentum diffusion tensor since the latter involve 
at most a second derivative of the correlation tensor [cf. 
Eqs.pl Ell pt)]. 



b. Interface contribution 

The interface-dependent part C\£ t of the field correla- 
tions is due to t wo c ontributions: the evanescent modes 
ff (k | , fi; r) (Eq.(Cl) with 1 < u < no), and the crossed 



term between waves incident from above and reflected at 



the interface (the two terms of fj_(k|,^;r) in Eq.(C8)). 
Collecting these contributions, one has 



C^(ri,r 2 ) = 

no 2tt 

10 M 
en(u,<p, ft)ey(u,tp, n)\t(u, fi)\ 2 x 
x exp[ik|| (r^i - r|| )2 ) + ik z (zi + z 2 )] 

1 2ir 



8- 



(eniu, V , n)e<$*(u, <p, tx)r*(u, /*) e -<M*i+*0 + 
+ eg>(«, <p, rie* u (u, ^ fx)r(u, fx) e <M*i+«d) x 
x exp[zk|| (r, Ul - r, U2 )] (C20) 

(Recall that lmk z > in the first integral and k z > 
in the second.) It is now evident that the interface con- 
tribution only depends on the in-plane difference vector 
s || = r ||.2 — r ||.i an d the sum of the distances Z\ + z%. 

We perform the integration over the azimuthal angle ip. 
with the help of the following formula |n] and its deriva- 
tives with respect to ks\\ 



dipe- lks * ucosv = 2irJ (k Sll u) 



(C21) 



where Jo is the Bessel function of zeroth order. We also 
observe the following property of the Fresnel coefficients 



n Q > u > 1 



u 2 -l 

2 9 



\t(u,n)\ 2 = 2Imr(u,/i) (C22) 



that allows to combine the contributions of reflected and 
evanescent modes in a compact way. Finally, the corre- 
lation tensor is decomposed into its isotropic, axial and 
quadrupolar parts, according to 

C^(z; S ||)=C^(r-I S ,r + Is) 

= Cg tfl (z; s,|) + C% t>1 (z; B „) + C^ )2 (z; 8 ,|) (C23a) 

C^ (z;s||)=c (z;sjj)^ (C23b) 

^ til (z;8||) = fca(z;sj|)(>'S - '\8 gd ) (C23c) 

C^ j2 (z;s||) = qo (z;sj)(s^S^ - + 

+fc 2 g 2 (z;s 2 )(sj |S j | - |sj}(5*J - <P' Z <F' Z )) (C23d) 

In these expressions, the dimensionless weight functions 
Co, qo, a, 1i are given by the Sommerfeld integrals 

1 f n ° du u 

c (z; s|) = - Re / J (ks\\u) x (C24a) 

2 J Q v 



x (txe + (2u 2 — 1)ttm) exp 2ikzv 



9o(^;s||) = - 



Re 



1 du u 
o v 



Jo{ksnu) X 



(C24b) 
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x (—i'te + {u 2 + 1)ttm) exp 2ikzv 

, 9N 3 f n ° , Mksuu) 
a(z;Sn) = -Im / duu — ; — - — r<rM exp 2i/cz?; (C24c) 
11 2 Jo 

q 2 (z;0)= S -Re ™ 
2 7o w 



(fc S ||)2 

x(rTE — (u 2 — 1)t*tm) exp 2ikzv 



(C24d) 



These integrals are computed numerically and are plotted 
as a function of fcz in Fig. [I] for sm = 0. The refractive 
index is no = 1.5. 



3. Relation to field susceptibility 

As a final comment, we would like to display the link 
between the vacuum correlation tensor on the one hand, 
and the classical field susceptibility, on the other. The 
latter quantity gives the (positive- frequency) electric field 
created at position r x by an oscillating dipole located at 
r 2 - 



S i (r 1 )e- iu! * = G^(r 1 ,r 2 )£» ie - 



(C25) 



The tensor G % J may be calculated from classical electro- 
dynamics p2|~|64|] . For the vacuum-dielectric interface, 
we have checked the following relation to the vacuum 
correlation tensor 



^(r 1 ,r 2 ) = 2/iImG ij '(r ll r 2 ). 



(C26) 



This result shows that, at least for zero temperature and 
a two-level system, both the spontaneous emission rates 
and the associated forces may be calculated from a clas- 
sical field calculation alone, witho ut explicitly quantizing 
the field. In particular, Eq.(C2£) justifies our interpre- 
tation of the radiation pressure force in terms of a clas- 
sical picture where the atomic dipole interacts with its 
own radiation reaction field reflected from the vacuum- 
dielectric interface. 



APPENDIX D: OPTICAL PUMPING IN A 
TM-POLARIZED EVANESCENT WAVE 

In this appendix, we study the optical pumping of a 
J = 1/2 atom in a n eva nescent wave with TM polariza- 
tion. The hclicity ( 87b ) is parallel to the y-axis that we 
choose as the quantization axis. The field is then a su- 
perposition of a + and er~ circular components, and the 
light-shift operator (|84|) is diagonal. It may also be veri- 
fied from the optical pumping equation Eq.(^6|) that the 
magnetization components J XiZ decouple from J y and the 
total population. Assuming that the atoms are initially 
unpolarized, the pumping process only depends on J y 
whose evolution is given by (we put again w = 1 in the 
optical pumping equation) 



8Jy 

dt 



2nQ 



T p (z) = 2a z T 



k 2 



T y ^ e i ax (z) = 2a 2 T' 00 e' 2KZ c ll (z) 



k 2 



(Dla) 
(Dlb) 
(Die) 



Optical pumping hence builds up a net magnetization 
along the negative y axis with a rate T p (z) that only de- 
pends on the fluorescence rate TooC^z) for a polarization 
parallel to the interface. This may be understood from 
the rate equations for the populations w± — i(l ± J y ) of 
the Zeema n sub levels | ± l/2) y with respect to the y axis. 
From Eqs.(Dl), the following rate equations are easily 
found 



dt 

-+(*) 



= -r_i — ,-(z)w + + r — , 
= — r — >+(z)w_ + r_i — ; 

2 .relax (^) -^p(^)] 
(«• 



2a 2 r^ e" 

.relax 

2a 2 T' e~ 



2 h. z 



.(z)w— 
-{z)w+ 

Q) 2 



2k 2 



' 2 -C||(z 



r,(2)] 

(/c-f 



Qf 



2k 2 



(D2a) 
(D2b) 

(D2c) 

(D2d) 



As expected, the rate for the transition | — 1/2) y — > | + 
l/2) y , e.g., is proportional to the square |£o+| 2 = \( K ~ 
Q)/\/2k\ 2 of the amplitude of the a + polarization in the 
evanescent field (the circular polarization being defined 
with respect to the y axis). In this process, the atom 
absorbs a er + polarized photon from the driving field and 
spontaneously emits a ir polarized photon (c/. Fig. f?]c). 
The electric field of the latter being parallel to the y 
axis, the spontaneous transition rate is proportional to 
the coefficient C y ' v (z) = cm(z) of the vacuum correlation 
tensor. 

The opp osite transition | + 1/2) y — ► | — l/2) y has a 
larger rate ( p2d| ) because in the TM case, the amplitude 
£o- = (n + Q)/\/2k of the o~ polarization is stronger 
than £o+- The steady-state magnetization is n ot m ax- 
imum, however, because the pumping rat e T v ( Dlb| ) is 



smaller than the relaxation rate T yre i ax ( Dlc| ), the po- 



larization of the TM evanescent wave being not purely 
circular. The steady-state solution of the Eq.(Dl) yields 



(stat) (stat) 
W — W±_ 



j(stat) 

-Jy 



^y. relax (-^) 

as stated in the text. 



2kQ 



< 1, 



(D3) 
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